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SUMMARY 

This thesis develops a computer solution to the exhaust gas 
reingestion problem for aircraft operating in the reverse thrust mode on a 
cross wind -free runway. The computer program determines the location 
of the inlet flow pattern, whether the exhaust efflux lies within the inlet 
flow pattern or not, and if so, the approximate time before the reversed 
flow reaches the engine inlet. The program is written so that the user 
is free to select discrete runway speeds or to study the entire aircraft 
deceleration process for both the farfield and cross -ingestion problems. 
While developed with STOL applications in mind , the solution is equally 
applicable to conventional designs . 

The inlet and reversed jet flow fields involved in the problem are 
assumed to be non-interacting. The nacelle model used in determining 
the inlet flow field is generated using an iterative solution to the Neuman 
Problem from potential flow theory while the reversed jet flow field is 
adapted using an empirical correlation from the literature. Sample 
results obtained using the program are included. 
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NOMENCLATURE 


Symbol 

Definition 


A' 1 

2 

Area , ft 


AL' 

Nacelle length, ft 


CHEKX 

Normal velocity on the end cap 


CHEKR 

Normal velocity on the nacelle 


D' 

Distance between two points , ft 


Dl' 

Distance between a point of interest in space and a 
on the engine inlet, ft 

point 

D2 ' 

Distance between a point of interest in space and a 
on the nacelle, ft 

point 

D3‘ 

Distance between a point of interest in space and a 
on the end cap, ft 

point 

DIAJET 

Diameter of the reversed jet at the point of origin 


DSAB 

Distance between points "a" and "b" 


EP 

Nacelle generation no-flow criteria - a select percentage 
of the freestream velocity 

g' 

Number of singularities per unit area, ft 


m' 

Strength of a singularity, ft^/sec 


PMPP 

The "P" coordinate in the jet plane of the Maximum 
Penetration Point of the reversed jet 



In this work, primes indicate dimensional quantities while non- 
primes indicate dimensionless quantities. Velocity and length terms are 
non-dimensionalized by referring them, respectively, to the freestream 
velocity U^' and to the nacelle radius R' . The product of U ' R’ is used 
to non-dimensionalize the velocity potential and stream function terms. 
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Symbol 

Q* 

QMPP 

qsl' 

r 

r3 

RA 

RAB 

RB 

RCROSS 

RMPP 

RPEST 

s j 

TAB 

U J 

V 

VAB 

VAX 

VAXY 


Definition 

3 

Volumetric flow rate from a singularity, ft /sec 

The "Q" coordinate in the jet plane of the Maximum 
Penetration Point of the reversed jet. 

3 

The strength of the inlet sink, ft /sec 

A radial space coordinate 

A radial space coordinate on the model nacelle 

Radius of point "a" 

Radial distance between points "a" and "b" 

Radius of point "b" 

Radius of the Maximum Penetration Point of the reversed 
jet of an adjacent engine on the coordinate system of 
an engine under study 

Radial coordinate of the Maximum Penetration Point of the 
reversed jet 

Radius of the pre-entry streamtube at an infinite axial 
location upstream of the engine inlet 

A velocity component of VELJET 

The time required for a tracer particle to go between points 
"a" and "b" f (TAB 1 ) (U^'j/CR') 

A velocity component associated with VELJET 

Velocity 

The average velocity of a tracer particle between points 
"a" and "b" 

A normal velocity produced by a singularity on the end cap 

The axial velocity at any point in space 
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Symbol 

VELJET 

VELRAT 

V ' 
n 

VP 

VPA 

VPB 

Vr 

VRAD 

w r 

X 

xl 

XA 

XAB 

XB 

XJET 

XMPP 

XOP 

XSPACE 

y 


Definition 

The velocity of the reversed jet 

The inlet-to-freestream velocity ratio ,qsl'/(2 • U ') 

A normal surface velocity, ft/sec 
Velocity at a point 
Velocity at point "a " 

Velocity at point "b" 

The normal velocity produced by a singularity on the 
nacelle 

The radial velocity of any point in space 

A velocity component of VELJET 

An axial space coordinate 

An axial space coordinate along the nacelle 

The axial coordinate of point "a" 

The axial distance between points "a" and "b" 

The axial coordinate of point "b" 

The axial coordinate of the reversed jet origin 

The axial coordinate of the Maximum Penetration Point 
of the reversed jet 

The axial coordinate of the Maximum Penetration Point of 
the reversed jet of an adjacent engine on the co- 
ordinate system of an engine under study 

The "x" component of the distance between the end cap 
center points of two adjacent engines 

A space coordinate 


x 



Symbol 


Definition 


yi 

YJET 

YMPP 

YOP 

YSPACE 

z 

zl 

ZJET 

ZMPP 

ZOP 

ZSPACE 

al 

a2 

P 

e . 

01 

'(> 


A space coordinate on the nacelle 

A coordinate of the reversed jet exhaust origin 

A coordinate of the Maximum Penetration Point of the 
reversed jet 

A coordinate of the Maximum Penetration Point of the 

reversed jet of an adjacent engine on the coordinate 
system of an engine under study. 

The "y" component of the distance between the end cap 
center points of two adjacent engines 

A space coordinate 

A space coordinate on the nacelle 

A coordinate of the reversed jet exhaust origin 

A coordinate of the Maximum Penetration Point of the 
reversed jet 

A coordinate of the Maximum Penetration Point of the 

reverS ed jet of an adjacent engine on the coordinate 
system of an engine under study 

The "z " component of the distance between the end cap 
center points of two adjacent engines 

The exhaust jet pitch angle 

The exhaust jet turning angle 

Angle between the reversed jet and the x-y plane 

The angle in the jet plane that the reversed jet makes with 
the "P" axis 

Angular space coordinate 

Velocity potential 

xi 



Symbol Definition 

col Velocity potential associated with the freestream-inlet 

combination 

^CAX Velocity potential associated with the distributed 

compensatory singularities on the end cap 

cp Velocity potential associated with the distributed 

compensatory singularities on the nacelle 


cp F g Velocity potential induced by the freestream 

cp Velocity potential induced by the inlet sink 

O 

Y Stream function 
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CHAPTER I 


INTRODUCTION 

The forthcoming development of both military and commercial 
short take-off and landing (STOL) aircraft will have a significant impact 
on air transportation. Quiet, civil STOL aircraft will greatly diminish 
city-to-city travel times and help relieve present airport congestion by 
operating from short, inner city airstrips. Such airstrips will be in- 
expensive enough that jet transportation can be extended to small 
communities and underdeveloped countries alike. In military versions, 
STOL transports will greatly improve the ability to supply remote regions . 

But before such a family of aircraft can be put into service there 
remain a number of problems to be resolved. One of the more important 
areas is the need for better thrust reversers. Unlike conventional air- 
craft, the normal means of stopping both military and commercial STOL’s 
will most likely be through the use of reversers alone. For the military, 
the capability to brake with reversers alone will greatly enhance 
operation from unprepared airstrips, by avoiding the rutting problem 
associated with wheel braking. For the commercial user, efficient 
reverse thrust braking is a matter of both economics and safety. The 
high operating costs of such aircraft will demand maximum daily 
utilization for profitable operation. Not only is brake maintainance a 
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major operating expense, but brake cooling requirements play a 
significant role in determining the aircraft turn-around time. From 
safety considerations, the full reverser stopping capability will improve 
operation from short, icy runways. 

Though thrust reversers have been used on jet aircraft for years , 
none have achieved the full stopping capability. There are several 
reasons for this. First, unlike STOL aircraft, conventional aircraft do 
not have sufficient thrust-to-weight ratios to stop within reasonable 
distances without the simultaneous use of wheel brakes. Secondly, the 
tendency of the engines to reingest their own reversed exhaust at low 
speeds has demanded that reverse thrust operation be terminated long 
before the aircraft has reached a halt. In addition to damaging parts , 
reingestion can cause compressor surge and greatly diminish the 
magnitude of the braking force. Thus any practical STOL aircraft must be 
designed to operate in reverse thrust, free of exhaust gas reingestion 
down to very low ground speeds. 

The purpose of this thesis is to develop an analytical model of 
the flow field near an aircraft engine operating in reverse thrust on a 
crosswind-free runway and to predict whether exhaust gas reingestion 
will occur. The study provides a tool, in the form of a computer program, 
to be used in the design of engine-nacelle-reverser systems. The 
designer is free to choose whatever dynamic and geometric conditions 
he wishes and can then see how effective they are with regard to the 
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reingestion problem. 

The overall flow field in this problem can be divided into two 
regimes: 

1 . the inlet flow field which is basically potential, 
and 

2. the reversed jet flow field which is highly 
turbulent. 

Ordinarily, the presence of two such diverse flow fields would make the 
development of any single analytical model an enormously difficult task. 
But in this investigation the two are mathematically uncoupled and 
solved separately. This greatly simplifies the analysis. 

There are several basic types of exhaust gas reingestion. The 
first type is near-field reingestion which occurs when the exhaust efflux 
passes too close to the nacelle. Because of the Coanda effect, the 
reversed jet attaches to the nacelle and subsequently enters the engine 
inlet. A second type is farfield reingestion where the jet penetrates the 
engine inlet flow pattern. A third type is cross -ingestion, where the 
reversed jet penetrates the inlet flow pattern of an adjacent engine. 

The computer program presented here is designed to study the 
latter two cases. It will not analyze the first case since no relavent 
studies of nacelle attachment currently exist. Therefore all solutions 
generated by this program are based on the assumption that the near- 
field problem does not occur. 
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The basic case under study is that of a turbofan engine installed 
in a long-duct nacelle, with a target reverser simultaneously handling 
the fan and core engine flows. One reason for selection of this con- 
figuration is its superiority in reducing approach and sideline engine 
noise. This currently is an important consideration since to receive 
community acceptance STOL aircraft will have to be significantly quieter 
than conventional aircraft. 



CHAPTER II 


PAST RESEARCH EFFORTS 

A search of the literature failed to reveal a past analytical 

solution to the thrust reverser reingestion problem. The solution 

presented here is based on a method proposed by Tatom [l ] . This 

method employs an axisymmetric model of an engine nacelle discharging 

2 

round, turbulent, reversed jets and is described in detail in the next 
chapter. 

Because of the importance of mathematically uncoupling the 
inlet and reversed-jet flow fields in Tatom's method, an investigation 
of the validity of this simplification was first made. It appears that 
the concept is well founded since there is ample evidence that: 

1 . The effect of the presence of the reversed jet 
on the freestream is small (i.e. , the freestream 
flow near the jet is essentially the same with 
and without the jet). 


1 Numbers in brackets indicate references cited in the 
Bibliography. 

2 

Jet effluxes from target type reversers tend to be approximately 
round in cross-section. 
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2. The effect of the engine inlet suction on the 
reversed jet is small (i.e., the trajectory of 
the jet is essentially the same with and without 
the presence of the inlet). 

Keffer and Baines [2] studied round turbulent jets introduced 
normally into a freestream. It was observed that the freestream was 
unaffected by the presence of the jet and that in the vicinity of the jet 
the static pressure and mean velocity of the freestream could be con- 
sidered constant. Additional evidence of non-interaction between the 
jet and the freestream can be found in the report by Weiss and 
McGuigan [3 ] . This paper contains oil-streak photographs of cold- 
flow tests of a model reverser-nacelle configuration. In each of these 
photographs the freestream is essentially parallel at a distance of less 
than 2 jet diameters upstream of the deflected jet. 

In each of the references cited above, the jets were discrete 
and thus produced relatively small blockage of the freestream flow field. 
The works of Cooper [4] and Hayden [5] are concerned with a two- 
dimensional reverser model with a simulated engine inlet. In this study 
the reversed jet blocked the entire freestream flow. Nowhere was the 
inlet flow field far removed from the reverser flow field . Yet flow 
visualization photographs of the reversed jet and temperature and velocity 
data taken with and without inlet suction, show no significant differences. 
These results further indicate that the freestream is almost unaffected 
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by the presence of the jet. Thus the independence of the two flow fields 
is verified. Since a round jet occupies much less volume in the region 
of the inlet than a two-dimensional jet, it follows that the effect of the 
inlet on a round jet will also be small. 

With the uncoupling hypothesis justified, attention is turned to 
previous efforts in the separate areas of inlet flow prediction and the 
trajectories of turbulent transverse jets. 

The inlet flow portion of the reingestion problem can be described 
adequately from potential flow theory. The basic problem (the Neuman 
Problem) involves generating a mathematical model of the flow near an 
engine nacelle within a freestream. This could be done by employing 
the Douglas -Neuman Potential Flow Computer Program developed by 
A.M.O. Smith and J. Pierce [6] . This selection was not used in [l ] 
for several reasons. First, the program was unavailable at the Vanderbilt 
University Computer Center and it was felt that adaptation of the program 
both to the Center's machine and to the reingestion problem would pre- 
sent as many difficulties as developing a new one. Secondly, it was 
felt that a new program might offer simplifications over the Douglas 
program and thus cost less to operate . 

The remaining area to be discussed is concerned with studies of 
jets penetrating into a freestream. Unfortunately, most of the available 
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literature is of little valve for two reasons: 

1. The studies are concerned with deflection of 
and velocities along the jet centerline which 
are of little importance here. The principle 
area of concern in the reingestion problem is the 
maximum penetration of the jet into the free- 
stream. 

2 . The sutdies are concerned with jet-to-freestream 
angles and velocity ratios considerably different 
than those encountered in reverser applications. 
Several analytical studies of opposing jets 
exist [7, 8, 9] , but these assume an 
incompressible, irrotational flow field. 

The engine-nacelle model proposed in [l ] incorporates the 
results of a Lockheed -Georgia study [10 ] to describe the reversed 
jets. The study was concerned with a round, turbulent jet introduced 
obliquely into an opposing freestream. Conducted in a low turbulence 
wind tunnel the experiment allowed the Maximum Penetration Point of the 
jet into the freestream to be photographically measured. The tests 
employed a wide range of values for the jet exit diameter, the jet-to- 
freestream velocity ratio, and the jet-to-freestream included angle. 

The result of this study was the Lockheed jet penetration correlation as 
presented in Figure 1, empirically relating the time averaged maximum 
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FIGURE 1. JET PENETRATION CORRELATION 
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penetration to the above menetioned variables . Additional data by 
Margeson [ll] is also presented . 

Figure 1 also shows a sketch of the reversed jet. Several 
characteristics of this jet should be noted. First, photographs indicate 
that the Maximum Penetration Point can be considered to lie approx- 
imately on an extension of the jet centerline. Secondly, the jet should 
not be considered as a fixed region in space. All flow visualization 
studies report that the reversed jet is an area of violent turbulence . 
Thus , the time-averaged Lockheed data does not show where the 
Maximum Penetration Point lies at any instant, but where it is most 


often found. 



CHAPTER III 


THEORY 

An overall view of the flow field involved in the problem is 
shown in Figure 2. As has already been noted, this flow field can be 
divided into inlet and reversed -jet flow fields; these being respectively, 
potential and turbulent in nature. The uncoupling hypothesis allows the 
two to be treated as independent problems. Hence the inlet flow model 
appears as a fictitious engine ingesting air but producing no exhaust, 
while the exhaust flow model appears as an isolated turbulent jet dis- 
charging obliquely into an opposing freestream. 

The inlet flow field is represented in the figure by the presence 
of the streamlines. Among these, the pre-entry streamtube is of special 
importance. This is defined such that fluid lying inside of it enters the 
engine while fluid lying outside of it travels past the nacelle. 

The exhaust flow too, has an item of special importance: the 
Maximum Penetration Point. At this point, the axial momentum of the jet 
has been completely depleted so that additional axial travel is determined 
by the opposing freestream. If the Maximum Penetration Point lies with- 
in the pre-entry streamtube, exhaust gas will be carried into the engine 
inlet. This is the cause of farfield reingestion. The primary task of this 
investigation is to develop an analytical model capable of generating 
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FIGURE 2. THE PHYSICAL SITUATION 
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the pre-entry streamtube and the Maximum Penetration Point, and 
determining whether or not the latter lies inside or outside of the 
former.* 

The Inlet Flow Field Model 

The inlet flow field solution is developed around a cylindrical 
nacelle configuration, with an end cap at the rear and a full frontal 
area inlet (Figure 2) . This geometry differs somewhat from actual 
engines . First, the inlet cannot fill the entire frontal area in a real 
nacelle due to structural and aerodynamic considerations. Secondly, 
nacelles are not cylindrical but are more streamlined bodies of 
revolution . The approximate nacelle representation is used because 
of the mathematical simplifications and resulting savings in computer 
time it affords. These simplifications are not believed to decrease 
significantly the accuracy of the solution since, in the farfield problem, 
the area where reingestion begins is somewhat removed from the engine. 
However, it must be recognized that near the inlet, the nacelle contour 
materially influences the shape of the streamlines. 


It should be noted that up to the Maximum Penetration Point the 
jet is entraining, not releasing, fluid. Hence, it is premissible for the 
reversed jet to lie within the pre-entry streamtube, as long as the 
Maximum Penetration Point does not. 
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Any exact mathematical model should acknowledge the rapid 
deceleration involved in the reversing process. While a transient 
description of the flow fields is desirable, a method for developing one 
is unclear. Therefore, as a final simplification, it is assumed that the 
aircraft passes continuously through a series of equilibrium flow fields 
in coming to rest. 

The purpose of the inlet flow field model is to generate the stream 
lines about the engine nacelle . The development of this model begins 
by placing a potential flow freestream (alligned with the axis) on an 
axisymmetric coordinate system to simulate the runway speed of the 
aircraft. The engine nacelle is generated within this freestream. 

Towards this end, an origin is established on the coordinate system 
and at this origin a disk sink is added to simulate the engine inlet 
(Figure 2) . 

In establishing the nacelle and end cap surfaces, a special set 
of boundary conditions must be satisfied. Because a real nacelle 
surface is solid, no fluid passes through it and hence the normal sur- 
face velocities must vanish. The same boundary condition applies along 
the end cap too, because the uncoupling assumption has removed the 
exhaust flow from the inlet model. The boundary conditions are satisfied 
by establishing a distributed system of compensatory singularities over 
the nacelle and end cap surfaces . 
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An iterative procedure is used to determine the singularity 

strengths. Initially, each singularity strength is set opposite and 

proportional to the normal velocity induced by the freestream-inlet 

combination at the point. This would be sufficient for compensation at 

isolated points . However, the presence of neighboring singularities 

induces an additional normal velocity at each point. These velocity 

components must also be cancelled and this is done by the iterative 

2 

adjustment scheme. Once the boundary conditions are satisfied the 

» 

inlet flow field model is capable of generating streamlines. 

The Reversed Tet Model 

The purpose of the reversed jet model is to locate the Maximum 
Penetration Point of the exhaust flow. The geometry of the problem is 
shown in Figure 3, with the nacelle outline included for clarity. The 
centerline of the jet is considered to lie in a plane. This jet plane (the 
P-Q coordinate system in the figure) is defined by the freestream and 
reversed jet velocity vectors and has its origin (o') at the point of 
reversed exhaust discharge. It is in this plane that the Lockheed 
correlation applies . 


2 

The surface generation process discussed here and proposed in 
[l ] evolved from the analysis of [ 1 2 ] . The major computational differ- 
ence between the two is that the latter makes no attempt at eliminating 
the normal velocity component induced by the neighboring singularities. 
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As the figure shows, the axisymmetric coordinate system of the 
inlet flow field model has been superimposed with a three-dimensional 
Cartesian system having the same origin. This new system is used to 
properly locate the origin of the jet plane with respect to the engine 
inlet. The angle between the freestream and the jet efflux in the P-Q 
coordinate system (0) is defined in terms of the pitch angle (cel) and 
turning angle (a2). The P-Q coordinates of the Maximum Penetration 
Point are found from the Lockheed correlation. Then, a multiple trans- 
formation of coordinates is employed to establish the axial and radial 
coordinates of the Maximum Penetration Point with respect to the 
original axisymmetric coordinate system. This completes the reversed 
jet flow field model. 

The geometries of the two flow fields are now superimposed to 
evaluate the likelihood of reingestion. If reingestion is predicted, the 
computer program is designed to note this and to determine the approx- 
imate time required for a fluid particle to travel from the Maximum 
Penetration Point to the engine inlet. 
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DEVELOPMENT OF THE MATHEMATICAL MODEL 

The Initial Singularity Strengths 
Before developing the inlet flow field model, it is useful to 
determine the initial compensatory singularity strength used in the 
iterative nacelle generation procedure . 

Consider an isolated area A' containing a singularity of strength 
m 1 . The volumetric flowrate Q' associated with this surface can be 
expressed as: 

Q' = m'A'g' 

where g' is the number of singularities per unit area and is equal to one 
in this case. In general, however, volumetric flowrate can be expressed 
as the product of a flow area and the velocity (V^ ) normal to it. For a 
singularity, the flow area is twice the area in the above equation because 
fluid simultaneously enters or leaves both sides. Hence: 

O' = 2A‘ • V' 
n 

Combining the two equations gives the singularity strength as a function 
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of V' , or: 
n 

m’ = 2 V‘ (for g'= 1) (1) 

n 

Therefore, the initial step in the iterative procedure is to set 
the strength of each singularity equal to twice the negative of the 
normal velocity induced by the freestream-inlet combination at the 
point to produce the cancelling normal velocity, V' . 

The Inlet Flow Field Model 

Let P in Figure 4 represent an arbitrary point in the flow field. 
The velocity potential induced at point P by all of the elements of the 
model can be described from potential flow theory as: 

- '"Vs + V S + ’’CR + '•’CAX (2a) 

where the terms on the right of the above expression represent the 
contributions to the potential from the freestream, the inlet sink, the 
distributed singularities on the nacelle, and the distributed singularities 
on the end cap, respectively. 

It is convenient to work with non-dimensional terms . Towards 
this end, length terms are non-dimens ionalized with respect to the 
nacelle radius R', while velocity terms are divided by the freestream 
velocity . The velocity potential terms are made dimensionless by 
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referring them to the product of the freestream velocity and nacelle 
radius . Equation 2a can thus be rewritten: 



U' R' 

CO 


*F8' + 


*S + 


cp 


CR 


+ cp 


CAX 


(2b) 


It is also convenient to group ther terms associated with the 
freestream and inlet sink together, or: 


cpl = 'Pps + 'Pg 
Hence, equation 2b becomes: 

cp = cpl + cp CR + cp CAX 


(3) 


(4) 


From potential flow theory, the radial and axial velocities at a 
point can be determined by taking the appropriate partial derivatives 
of equation 4; e.g. 


and 


VRAD = 


VAXY = 


d cp 

d cpl 

^CR 

3CP CAX 

~ dr 

dr 

+ 8r + 

dr 

dcp 

dcpl 

Stp CR 

8r|, CAX 

dx 

dx 

+ dx 

dx 


(5) 

(6) 
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The four velocity potential terms in equation 2b are now to be 
developed. From potential flow theory, a freestream can be described 
as: 

cp’s = ir x' 

or in non-dimensional terms: 

®FS =X (7) 

The remaining terms in equation 2b describe surfaces of 
distributed singularities. Figure 5 shows an arbitrary surface divided 
into subareas dA 1 , each containing one singularity of strength m'. 

The incremental velocity potential induced at a point P by any such 
subarea, a distance D' away can be described by [ 1 3 D : 


dcp ' 


1 m'dA' 
4 tt D' 


Hence, the velocity potential induced at P by all such subareas can 
be written as: 


<P 


- f 

\ TT J . 


m'dA' 


4n J D' 
A 


or in non-dimensional form: 


1 r m'dA 1 

4TT U' R' A D' 
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If, in Figure 5, point P lies on the surface, then the subarea con- 
taining P must be excluded from the integration of equation 8 . 

Two of the components of the model to be described by equation 
8 are disks. Hence as Figure 4 shows the area integral in equation 8 is 
evaluated with respect to r3' and Ol and can be written: 

dA' = r3' dr3' del 
or 

dA' = (R') 2 r3 dr3 d61 (9) 

The remaining component to be described by equation 8 is the cylinder 
and hence the area integral is evaluated with respect to xl l and 01; or: 

dA’ = R' dxl' d91 
or 

dA' = (R') 2 dxl d 01 (10) 

With continuing reference to equation 8, Figure 4 shows that 
there are three D' terms: one for the inlet (Dl'), one for the cylinder 
(D2‘) and one for the end cap (D3') . In general, the distance between 
point P and any point on the nacelle can be written: 

D' =\/ (x'-xl 1 ) 2 + (y'-yl') 2 + (z 1 - zl') 2 
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or 

D' = (R')\/(x-xl) 2 + (y-yl)^ + (z - zl) 2 (11) 

The terms x, y, and z, are the coordinates of point P from the 
origin. The terms xl , yl , and zl , are the coordinates of any point on 
the nacelle. Because the system under study is axisymmetric , reference 
point P can be defined as always lying at z = 0. 

In evaluating the Dl' and D3‘ expressions, the areas involved 
are disks and hence the xl terms are zero for the former and AL for the 
latter. In all the D' expressions, the values of yl and zl can be 
written: 

zl = r3 sin 0 

and 

yl = r3 cos 0 

In the Dl' and D3' expressions, r3 terms in the above pair of 
equations remain as variables while in the D2' expression r3 is a 
constant with a value of unity. The three distance equations can thus 
be written, after simplification, as: 

Dl’ = (R')\/x 2 + r 2 - 2r r3 cos 01 + r3 2 (12) 

D2' = (R')\/ (x-xl) 2 + r 2 - 2r cos 01 + 1 (13) 

D3‘ = (R')\/ (x-AL) 2 + r 2 - 2r r3 cos 01 + r3 2 


and 


(14) 
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The remaining term to be considered in equation 8 is the 
singularity strength, m'. For the case of the inlet sink, the strength 
(qsl') is constant over the area and is a function of a particular engine 
design and/or engine power setting. Thus it can be brought outside the 
integral . 

In the cases of the nacelle and end cap surfaces, m' is a function 
of axial position along the nacelle and radial position along the end 
cap. Therefore in these cases m* must remain inside the integrals. As 
mentioned in the previous section, the strength of each singularity on 
these surfaces is initially set at twice the negative of the normal 
velocity induced by the frees tream -inlet combination at the 
point. 

All of the expressions necessary to describe the velocity poten- 
tial at an arbitrary point can now be written. Combining equations 7, 
8,9, and 12 with equation 3 gives, after simplification, the potential 
due to the freestream and inlet: 


2l7 1 


r3 dr3 d 81 

[x 2 + r 2 - 2r r3 cos 01 + r3 2 ] * 


where the limits of integration in this and all of the velocity potential 
expressions are those already noted for the surface areas under con- 
sideration. 

An expression for the potential due to the nacelle can be obtained 
by combining equations 1, 8, 10, and 13 and gives, after simplification: 
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cp, 


CR 


ZTT 

_ j_ f r 
2tt J J 


2TT ~ AI Vr(xl) dxl del 

2 0 r ~" 1/2 
o o [(x-xl) + r - 2r cos 81 + 1 ] 


(16) 


And finally, combining equations 1, 8, 9, and 14 gives, upon 
simplification, the potential due to the end cap: 


2n .1 


<T> 


CAX 


J 


VAX (r3) r3 dr3 d01 


o o 


9 9 2 1 /2 

[(x-AL) + r - 2r r3 cos 01 + r3 ] 


(17) 


In order to evaluate equations 5 and 6 the partial derivatives of 
the above three equations must be taken with respect to x and r. 
Applying Leibnitz's rule [14 ] to equations 15, 16, and 17 gives, upon 
simplification: 


2tt 1 


(r - r3 cos 91) r3 dr3 d 61 


Bcpl _ qsl r r 

dr 4 tt J J r 2 2 „ . .. ^ -2,3/2 

o o [x + r - 2r r3 cos 01 + r3 J 


(18) 


dr tR 1 f‘ 2TT f AL [Vr (xl)][r - cos 61 ] dxl dOl 

8r 2 tt J J r 2 2 _ Q1 , 3/2 

o o [(x-xl) + r - 2r cos 81 + 1 J 


(19) 


drp 


2tt _1 


"CAX J_ r~" P* [VAX (r3)][r - r3 cos 91 ] [r3 ] dr3 d 61 

r 2Tr o o [(x-AL) 2 + r 2 - 2r r3 cos 01 + r3 2 ] 3//2 


( 20 ) 


dcpl _ qsl .x 


2 TT 1 


Bx 1 4 


j 


r3 dr3 d 6 1 


o v o [x 2 + r 2 - 2r r3 cos 8 1+ r3 2 ] 3/ ^ 2 


( 21 ) 
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Sep 

3x 


CR 


2tt 


2rr AL 

J J 

o o 


[Vr (xl) ][x-xl ] dxl d0 1 

5 o — T7T 

[(x-xl) + r - 2r cos 61 + 1 ] 


( 22 ) 


d^CAX 

Sx 


J_ f 2TT f 1 [VAX (r3)]fc-AL]&-3] dr3 d6l 

2TT "o o [(x-AL) 2 + r 2 - 2r r3 cos 01 + r3 2 ] 3/2 


(23) 


The above equations must be integrated to obtain the six 
velocity components . Integrations involving both variables of equations 
19, 20, 22, and 23 must be performed numerically because of the 
dependence of the singularity strengths on position. Equations 18 and 
21 , however, can be integrated in closed form with respect to r3, though 
they must be integrated numerically with respect to 01 . Performing the 
cosed form integration gives: 


,2tt 


Scpl _ qsl r \F a „r (2B + 4A) ~| 

ar - - TF J Iv-r-aJL 4 ^ - -== J 


4A-B 


V A + B + 1 


- cos 81 


{( 


_1)( 2B -4A + 2AB . 2bA \ 
4A-B^ 'V A + B + 1 


+ i 


■[* 


* + d01 
JR + B/2 “* 


(24) 
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and 




(2B + 4A) 


4A+B 


VA + B + 1 


•J 


2 2 

where A = x + r 


and 


B = -2r cos 01 


d01 (25) 

(26) 

(27) 


Numerical results using equations 5 and 6 can be obtained when 
the singularity strengths have been evaluated. 


Evaluation of the Singularity Strengths - 
The Iterative Procedure 


To satisfy the no-flow condition the normal velocity must vanish 
at each nacelle and end cap singularity. As explained earlier, the first 
step toward this end is to set the singularity strength at each location 
equal to twice the negative of the normal velocity induced by the free- 
stream-inlet combination. Consider any point P on the nacelle surface, 
The normal velocity CHEKF at P is determined from the equation: 


CHEKR = 


dcpl 

dr 


dcp 


- Vr + 


CR 


dcp 


CAX 


dr 


dr 


(28) 


at P 


In the above equation, the quantity Vr represents the velocity 

produced by the singularity at P and is initially equal in magnitude to 

. The remaining terms are identical to the corresponding terms in 
dr 

equation 5. However, because point P lies on the nacelle surface, the 
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area containing P must be excluded from the numerical integration of 

the ^ c °CR term (equation 19) . 

Tr 

The no-flow boundary condition is satisfied along the nacelle 
when the velocity CHEKR. vanishes at each singular point. 

An equation analogous to equation 28 for the normal velocity at 


any point P on the end cap is: 


dcpl 

CHEKX = VAX + 

ax 


Scp CR 8lp CAX 
&x + dx 


at P 


(29) 


The partial differential terms in the above equation are identical 

to the corresponding terms in equation 6. The quantity VAX represents 

the velocity produced by the singularity at P and is initially equal in 

magnitude to . In performing the numerical integration of the 

a x 

d fp CAX term, the area containing P must be excluded, 

ax 


The no-flow boundary condition is satisfied along the end cap 
when the velocity CHEKX vanishes at each singular point. 

The iterative procedure used in generating the nacelle and end 
cap surfaces is as follows: 

1. The normal velocity induced by the freestream- 
inlet combination is calculated at each nacelle 

("Ir"^) anC * enC * ° ap ("lx”) s ^- ngu ^ ar P°i- nt by 
evaluating equations 24 and 25, respectively. 
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2. At each nacelle singular point, Vr is set 
equal to the negative of , while at 
each end cap singular point, VAX is set 

, acpi 

equal to the negative of — — . 

o X 

3 . The normal velocity CHEKR is determined at 
each nacelle singular point by evaluating 
equation 28 and the normal velocity CHEKX is 
determined at each end cap singular point by 
evaluating equation 29. 

4. At each nacelle singular point where CHEKR 
is non-zero, an adjustment scheme (to be 
described below) resets Vr. Likewise, at 
each end cap singular point where CHEKX is 
non-zero, VAX is reset. 

5 . The procedure is repeated from STEP 3 until the 
ho-flow condition is met at all singular points 
to some specified precision. 

The adjustment procedure mentioned in STEP 4 is as follows. 
Consider first equation 28, Figure 6 represents all of the velocity 
components associated with this equation for a singular point P on 
the nacelle surface. Let it be assumed that CHEKR at P is not zero 
and hence Vr must be adjusted. The first step is to determine if 




VELOCITY COMPONENTS AT POINT P 


POINT P 



FIGURE 6. SINGULAR POINT VELOCITY VECTORS 


GO 

N) 
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CHEKR and Vr are of the same sign. If they are not, the magnitude 
of Vr must be increased. This is done by employing the equation: 

Vr 1 = Vr j ,, - CHEKR /2 (30) 

[new fold 

The choice of the correction term, CHEKR/2, in equation 30 is arbitrary. 

If Vr and CHEKR are of opposite sign, the magnitude of Vr 
must be decreased. The procedure for doing this depends on the value 
of CHEKR/2. If the magnitude of CHEKR/2 is less than that of Vr , 
equation 30 is applied. However, if the magntide of CHEKR/2 is greater 
than that of Vr, the new Vr is obtained from: 

Vr I = Vrl / 2.00 (31) 

'new 'old / 

Again, the choice of the correction is arbitrary. 

An identical adjustment procedure is used along the end cap with 
CHEKX (from equation 29) substituted for CHEKR and VAX substituted for 
Vr in equations 30 and 31. 

Generation of the Streamlines 

Let a fluid particle be released from a point P in the flow field 
with the object being to determine the path it follows . Since the flow 
field is assumed to be in equilibrium, the fluid particle will travel 
along a streamline. The value of the stream function along any 



34 


streamline is constant and hence: 


3Y 



streamline 


(32) 


In an axisymmetric system, the stream function is a function of 
x and r and hence: 

Y = Y (x,r) 

Taking the derivative and applying it along a streamline gives: 


3Y 



dr + 


3Y 

ax 


dx 


or rearranging: 


Y = Const. 


dr 

dx 


aY /dx 

aY/a r 


(33) 


The radial and axial velocity at any point in an axisymmetric 
system can be expressed from Stokes Stream Function [15] as: 


VAXY = - 

r 3r 


and 

VRAD = 


I 

r ax 


Rearranging the above two equations and combining them with 
equation 33 gives, after simplification: 


dr 

dx 


Y = Const. 


VRAD 

VAXY 


( 34 ) 
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With the coordinates of the release point known, the above 
differential equation can be numerically solved using the Runge- 
Kutta Method [ 1 6 ] to determine the radial and axial coordinates of 
the points along the streamline which the fluid particle follows . The 
velocities VRAD and VAXY needed in this solution are obtained by 
evaluating equations 5 and 6, respectively. 

The Reversed let Flow Field Model 

This section developes the axial and radial coordinates of the 
Maximum Penetration Point of the reversed jet relative to the 
axisymmetric coordinate system of the inlet flow field model. The 
problem is defined in Figure 3. 

The P coordinate of the Maximum Penetration Point in the jet 
plane is found from the Lockheed correlation (Figure 1): 

PMPP = 2.97 (DIAJET) (VELJET) ‘ 94 (1 - . 734SIN ’ 685 0) (35) 

From Figure 3, the Q coordinate of the Maximum Penetration Point can 
be written: 

QMPP = PMPP tan 0 (36) 

The axial coordinate of the Maximum Penetration Point with 
respect to the axisymmetric coordinate system can be written from 
Figure 3 as: 

XMPP = XJET - PMPP 


(37) 
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where PMPP is found from equation 35. An expression must now be 
developed for the radial coordinate of the Maximum Penetraion Point 
with respect to the axisymmetric coordinate system. From Figure 3: 

RMPP A/ YMPP 2 + ZMPP 2 (38) 

where; 

YMPP = YJET + QMPP • cos ((3) 

and 

ZMPP = ZJET + QMPP • sin (3) 


Substituting the above expressions into equation 38 gives, 
upon simplification: 


RMPP ==\/ YJET 2 +ZJET 2 + 2 QMPP [YJET-cosp +ZJET -sin(3 ] + QMPP 2 

(39) 

Everything necessary to solve equation 39 has now been 
developed except for the in-plane angles, Band 3. From Figure 3: 


or: 


Sj = VELJET • cos al 


J 


VELJET 
Likewise , for u 


= cos al 


J 


or: 


Uj = Sj cos a2 


U J 

— = cos a2 

S J 


(40) 


( 41 ) 
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From Figure 3: 

cos 9 = Uj AELJET H2) 

Equation 42 can be rearranged as: 

S J U J 

COS 9 " VELJET ' Sj 

Combining equations 40 and 41 with the above gives: 
cos 0= cos ocl • cosa2 


or; 


and 


0 = cos -1 [cos otl . cos a.2 ] 

Referring to Figure 3 , the following can be written: 
tan al = w j/ s j 

sin a2 = v j/ s j 


Also , 


tan P = Wj/Vj 


Equation 46 can be rearranged as: 


w s 

tan P = — • — 

S I. V J 


(43) 


(44) 

(45) 


( 46 ) 
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Combining equations 44 and 45 with the above gives: 


tan (3 = tan al/sin ct2 


or: 


-1 r TAN al n 

P = TAN — -T7C ! 

L sin cl2 


(47) 


With equations 43 and 47 complete, the location of the Maximum 
Penetration Point with respect to the axisymmetric coordinate system 
can be determined from equations 3 7 and 39. 


Time Calculations 


This section develops a method of approximating the time required 
for a fluid particle to travel from the Maximum Penetration Point to the 
engine inlet, for cases where reingestion occurs. Figure 3 shows two 
points, A and B, on a streamline located within the pre-entry streamtube. 
Equations have already been developed for finding both the radial and 
axial velocities at these points (equations 5 and 6), as well as the 


locations of the points themselves, (equation 34). Let the quantity VAB 
be defined as the average speed between points A and B. It follows 


then, that: 


VAB = 


VPA + VPB 
2 


(48) 


where VPA and VPB are the speeds at points A and B, respectively. 
These speeds can be determined from: 


VP A/ (VAXY) 2 + (VRAD) 2 


(49) 
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Let XAB be the axial distance between the points and RAB be the radial 
distance between the points. Then: 


XAB = XA - XB 


and 


RAB = RA - RB 


(50a) 
(5 0b) 


The approximate distance between the points, DSAB can be described as 


DSAB 


A/ (XAB) 2 + (RAB) 2 


(51) 


From the elementary equation, distance equals speed times time, 
DSAB can also be expressed as: 

DSAB = (VAB) (TAB) 


where TAB is the time required for the particle to travel the distance 
DSAB at an average speed of VAB. The above equation can be rewritten 


as; 


TAB = 


DSAB 

VAB 


(52) 


By summing the TAB values between all of the points along the stream- 
line from the Maximum Penetration Point to the inlet, the approximate 
time involved in the reingestion process can be determined. 


Cross Ingestion 

Figure 7 shows a sketch of a four-engined jet transport with 
wing-mounted engines . The quantities XOP and RCROSS are defined 
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respectively as the axial and radial coordinates of the Maximum 
Penetration Point of the reversed jet of the inboard engine relative to 
the outboard engine coordinate system. Clearly cross ingestion occurs 
if RCROSS lies within the outboard engine pre-entry streamtube at 
XOP. With this in mind, expressions for these quantities are now 
developed. 

The first step is to establish the center point of the inboard 
engine end cap relative to the same point on the outboard engine with 
the quantities XSPACE, YSPACE, and ZSPACE. With this done, XOP can 
be described as: 

XOP = XMPP - XSPACE (53) 

From Figure 7, the following expressions can also be written: 

YOP =• YSPACE - YMPP (54) 

and 

ZOP = ZSPACE + ZMPP (55) 

From these two equations, RCROSS can be described as: 

RCROSS A/yOP 2 + ZOP 2 (56) 

With XOP and RCROSS determined, the likelihood of cross 
ingestion can be evaluated. 
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The Pre-entry Streamtube Radius Equation 
The most important streamtube is the pre-entry streamtube and 
while the radius of this boundary cannot be analytically determined at 
arbitrary axial locations, it can be determined at minus infinity (RPEST). 
At minus infinity, the velocity within the pre-entry streamtube is the 

,K \ 

freestream velocity = 1 y . At the engine inlet, the velocity 

CO 

within the pre-entry streamtube is the sum of the freestream and inlet 

induced (VELRAT) velocities . The flow areas at these two axial locations 
2 2 

are rr . RPEST and rr . 1 , respectively. Since continuity is main- 
tained within the pre-entry streamtube, the following can be written: 


V • A = V . A 

at - 00 at Inlet 


2 2 
1 • TTRPEST = (VELRAT + 1) TT . 1 


or: 


RPEST 


A/ 


VELRAT + 1 


( 57 ) 



CHAPTER V 


THE COMPUTER PROGRAM 

The statement listing of the computer program is presented in 
Appendix I. The block diagram of this program is shown in Figure 8 and 
for clarity, the step numbers in the figure are included in the statement 
listing. 

Referring to Figure 8, the initial step in the program is the 
inputting of the geometric, dynamic, and program variables and selection 
of the program options . The procedure for doing this is described in 
Appendix II. 

Program initialization for the first inlet-to-freestream velocity 
ratio to be studied occurs in STEP 2 . 

In STEP 3 , the normal velocities induced at each nacelle and end 
cap singularity by the freestream-inlet combination are computed from 
equations 24 and 25, respectively. 1 In addition, the initial strength 
settings of the compensatory singularities are assigned here. 

The iterative nacelle generation procedure comprises STEPS 4,5, 
and 6. In STEP 4, the induced normal velocities at each nacelle and end 

■^For programming convenience, the number of singularities on the 
nacelle and end cap are equal. 
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acpi 
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Scpl 
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INITIAL STRENGTH SETTINGS OF 
COMPENSATORY SINGULARITIES 


dM 


i CALCULATION OF CHEKR 


n 
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NO 
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\ 
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\ 

WHAT DO THE CALCULATIONS SHOW? 
^ 


@- j HAVE ALL VELOCITY RATIOS BEEN STUDIED? 
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1 
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FIGURE 8 . FLOW CHART OF THE COMPUTER PROGRAM 
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cap singularity due to all of the elements of the system are computed 
from equations 28 and 29, respectively. The purpose of STEP 5 is to 
record in computer memory those points where the above normal 
velocities have not vanished. At all such points the singularity 
strengths are adjusted in STEP 6 using the procedure described in 
Chapter IV. The degree of accuracy to which the no-flow condition is 
established is controlled by the inputted quantity EP, which represents 
a selected percentage of the freestream velocity. 

The program proceeds to STEP 7 when the no-flow condition is 
met at every singular point. Otherwise, control is returned to STEP 4. 

The purpose of STEP 7 is to provide the coordinates of a starting 
point for the streamline generation procedure of STEP 8 . Two program 
options are available here. With one option the coordinates of the 
Maximum Penetration Point of the reversed jet are computed using the 
Lockheed correlation. This option is employed to evaluate the likelihood 
of reingestion. If the other option is chosen, the coordinates of an 
inputted point are used. This option is employed for generating 
selected streamlines . 

With an initial point determined, the path of a streamline is 
generated in STEP 8. Additionally, the time required for a fluid particle 
to travel the streamline is determined here. 

In STEP 9 the likelihood of reingestion is evaluated. Output 
confirms whether or not exhaust efflux has entered the engine inlet. If 
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reingestion occurs, the fluid particle time is also outputted. Addition- 
ally, a program option is available to determine if the entraining portion 
of the jet penetrates the pre-entry streamtube. 

Program operation is terminated in STEP 10 unless further inlet- 
to-freestream velocity ratios are to be studied. If this is the case, 
control is transferred back to STEP 2. 



CHAPTER VI 


RESULTS AND DISCUSSION 

Varification of the operationally of the computer program con- 
sists of three steps: 

1. Demonstration of the nacelle generation sections 
of the program. This is accomplished by showing 
that the inlet flow field model can be generated 
over a wide range of geometric and dynamic 
conditions, to any specified degree of accuracy. 

2. Demonstration of the streamline computational 
scheme. This is accomplished qualitatively by 
plotting selected streamlines about a nacelle and 
quantitatively by showing that continuity is 
satisfied between adjacent streamlines. 

3. Demonstration of the ability of the program to 
analyze a realistic reingestion problem. 

Nacelle Generation 

Results from several nacelle generation studies are presented in 
Tables 1 through 5. At each singularity, the tables note the normal 
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TABLE 1 . NACELLE GENERATION DATA 


VELRAT = 60; ASPECT RATIO = 3; 

EP = 0.01; 

NO. OF ITERATIONS = 12 


Singularity 

T tinn 

dcpl 

dr 

Vr 

CHEKR 

dcpl 

dx 

VAX 

CHEKX 

i 

-29.12381 

24.17546 

0.007 

0.18362 

-0.75719 

0.010 

2 

-25.46819 

19.73041 

0.005 

.18371 

- .75737 

0.010 

3 

-14.12882 

9.91636 

0.003 

.18396 

- .75785 

0.010 

4 

- 8.70776 

4.58291 

0.001 

.18436 

- .76371 

0.005 

5 

- 5.66953 

2.44176 

-0.010 

.18494 

- .76497 

0.005 

6 

- 3.84147 

1.35647 

-0.010 

.18566 

- .76272 

0.010 

7 

- 2.69018 

.78559 

-0.009 

.18656 

- .76994 

0.005 

8 

- 1.93872 

.47825 

-0.007 

.18760 

- .77259 

0.006 

9 

- 1.43300 

.30737 

-0.005 

.18879 

- .77628 

0.007 

10 

- 1.08323 

.20862 

-0.004 

.19016 

- .78017 

0.008 

11 

- 0.83533 

.14861 

-0.004 

.19168 

- .78562 

0.009 

12 

- .65568 

.10929 

-0.005 

.19334 

- .79736 

0.005 

13 

- .52286 

.09744 

0.008 

.19518 

- .80610 

0.006 

14 

- .42287 

.08309 

0.009 

.19716 

- .81636 

0.008 

15 

- .34636 

.07186 

0.006 

.19929 

- .83018 

0.010 

16 

- .28693 

.06938 

0.005 

.20154 

- .85308 

0.007 

17 

- .24015 

.07722 

0.005 

.20397 

- .87755 

0.008 

18 

- .20286 

.08830 

-0.008 

.20653 

- .91593 

0 .005 

1 9 

- .17282 

.14561 

-0.007 

.20923 

- .96087 

U . 008 

20 

- .14835 

.29073 

-0.008 

.21207 

-1.01060 

0.010 

21 

- .13781 

.50227 

-0.007 

.21354 

-1.03314 

0.005 


4 ^. 

oo 



TABLE 2 . NACELLE GENERATION DATA 



VELRAT 

= 10; 

ASPECT RATIO = 3; 

EP = 0.01; 

NO. OF ITERATIONS = 12 


Singularity- 


ScdI 

Vr 

CHEKR 

Seel 

VAX 

CHEKX 

Location 


Sr 



Sx 



1 

-4 

.85397 

4.03431 

0.007 

0.86394 

-1.01892 

0.007 

2 

-4 

.24470 

3.29298 

0.006 

.86395 

-1.01914 

0.007 

3 

-2 

.35480 

1.50842 

0.007 

.86399 

-1.01982 

0.007 

4 

-1 

.45129 

0.76922 

0.006 

.86406 

-1.02095 

0.007 

5 

-0 

.94492 

.40662 

-0.003 

.86416 

-1 .02260 

0.007 

6 


.64025 

.22387 

-0.006 

.86428 

-1.02412 

0.008 

7 

_ 

.44836 

.12865 

-0.006 

.86443 

-1.02690 

0.009 

8 

_ 

.32312 

.07572 

-0.008 

.86460 

-1.03038 

0.009 

9 

— 

.23883 

.04957 

-0.006 

.86480 

-1.03468 

0.010 

10 


.18054 

.03204 

-0.008 

.86503 

-1.04505 

0.006 

11 

_ 

.13922 

.02505 

-0.005 

.86528 

-1.05205 

0.006 

12 

— 

.10928 

.01624 

-0.009 

.86556 

-1.05957 

0.008 

13 

— 

.08714 

.01584 

-0.007 

.86586 

-1.07019 

0.009 

14 

■■ - 

.07048 

.01789 

-0.004 

.86619 

-1.08857 

0.005 

15 

— 

.05773 

.02283 

-0.002 

.86655 

-1.10572 

0.006 

16 

_ 

.04782 

.03210 

-0.001 

.86692 

-1.12854 

0.008 

17 

_ 

.04002 

.04972 

-0.002 

.86733 

-1.15886 

0.009 

18 


.03381 

.07909 

-0.010 

.86776 

-1.20559 

0.007 

19 

_ 

.02880 

.16141 

-0.006 

.86820 

-1.26252 

0.009 

20 

_ 

.02472 

.35474 

-0.009 

.86868 

-1.33049 

0.005 

21 

- 

.02297 

.63382 

-0.008 

.86892 

-1.35254 

0.006 



TABLE 3. NACELLE GENERATION DATA 



VELRAT = 5; ASPECT RATIO = 2; 

EP = 0.01; 

NO. OF ITERATIONS = 14 


Singularity- 

acoi 

Vr 

CHEKR 

dcpl 

VAX 

CHEKX 

Location 

dr 



dx 



1 

-2.42699 

2.00175 

0.009 

.85071 

-1.01983 

0.008 

2 

-2.74531 

2.25474 

0.009 

.85075 

-1.02006 

0.008 

3 

-1.71072 

1.22803 

0.007 

.85084 

=1.02075 

0.008 

4 

-1.17740 

0.74237 

0.004 

.85099 

-1.02192 

0.008 

5 

-0.84646 

.47211 

0.005 

.85121 

-1.02360 

0.009 

6 

- .62576 

.30499 

0.002 

.85150 

-1.02485 

0.010 

7 

- .47246 

.20102 

-0.001 

.85184 

-1.03277 

0.005 

8 

- .36304 

.13840 

0.001 

.85224 

-1.03664 

0.006 

9 

- .28331 

.09497 

-0.002 

.85271 

-1.04080 

0.007 

10 

- .22418 

.06794 

-0.004 

.85323 

-1.04632 

0.008 

11 

- .17967 

.05194 

-0.005 

.85381 

-1.05296 

0.009 

12 

- .14569 

.04375 

-0.005 

.85445 

-1.06194 

0.010 

13 

- .11942 

.04371 

-0.002 

.85514 

-1.07872 

0.006 

14 

- .09886 

.04334 

-0.006 

.85589 

-1.09307 

0.008 

15 

- .08260 

.05370 

-0.006 

.85669 

-1.11731 

0.005 

16 

- .06961 

.06959 

-0.010 

.85754 

-1.14520 

0.007 

17 

- .05912 

.10802 

-0.005 

.85844 

-1.18936 

0.005 

18 

- .05058 

.16379 

-0.008 

.85939 

-1.25380 

0.008 

19 

- .04357 

.27233 

-0.006 

.86039 

-1.37106 

0.007 

20 

- .03777 

.51670 

-0.009 

.86143 

-1 .56440 

0.005 

21 

- .03524 

.91305 

-0.010 

.96197 

-1.66026 

0.007 



TABLE 4. NACELLE GENERATION DATA 



VELRAT = 5; ASPECT RATIO = 4; 

EP = 0.01; 

NO. OF ITERATIONS = 11 


Singularity- 

Sep 1 

VR 

CHEKR 

6 CD 1 

VAX 

CHEKX 

Location 

3r 



ox 



1 

-2.42699 

2.02620 

0.008 

.96139 

-1 .05643 

0.007 

2 

-1.71072 

1.25887 

0.007 

.96139 

-1.05665 

0.007 

3 

- .84646 

.48609 

0.006 

.96140 

-1.05732 

0.007 

4 

- .47246 

.21190 

0.005 

.96141 

-1.05844 

0.007 

5 

- .28331 

.09349 

-0.004 

.96142 

-1.06006 

0.007 

6 

- .17967 

.04307 

-0.008 

.96144 

-1.06222 

0.008 

7 

- .11942 

.02236 

-0.007 

.96147 

-1.06497 

0.008 

8 

- .08260 

.01461 

-0.004 

.96150 

-1.06841 

0.008 

9 

- .05912 

.00929 

-0.004 

.96153 

-1.07265 

0.009 

10 

- .04357 

.00198 

-0.009 

.96156 

-1 .07781 

0.010 

11 

- .03293 

.00170 

-0.008 

.96161 

-1.08912 

0.005 

12 

- .02543 

.00245 

-0.005 

.96165 

-1.09688 

0.006 

13 

- .02001 

.00400 

-0.003 

.96170 

-1.10670 

0.007 

14 

- .01601 

.00644 

0.000 

.96176 

-1.11870 

0.008 

15 

- .01300 

.01300 

0.006 

.96182 

-1.13342 

0.009 

16 

- .01069 

.01624 

0.005 

.96188 

-1.15675 

0.005 

17 

- .00889 

.02644 

0.005 

.96195 

-1.17934 

0.007 

18 

- .00747 

.04486 

-0.002 

.96202 

-1.20706 

0.008 

19 

- .00633 

.10096 

-0.006 

.96209 

-1.23760 

0.009 

20 

- .00541 

.26865 

-0.009 

.96217 

-1.26875 

0.005 

21 

- .00502 

.51886 

-0.008 

.96221 

-1.27575 

0.006 
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TABLE 5 . EFFECTS OF NO-FLOW CRITERIA 


VELRAT = 5; ASPECT RATIO = 3 



PART A 
EP =0.01 

NO. OF ITERATIONS = 12 

PART B 
EP = 0.001 

NO . OF ITERATIONS = 1 9 

Singularity 

Location 

Scpl 

Br 

CHEKR 

Vr 

Vr 

CHEKR 

1 

-2.42699 

0.008 

2.02155 

2.01511 

0.001 

3 

-1.17740 

0.004 

0.75460 

0.75191 

0.000 

5 

-0.47246 

-0.002 

0.20368 

0.20418 

-0.001 

7 

-0.22418 

-0.007 

0.06212 

0.06643 

-0.001 

9 

-0.11942 

-0.006 

0.02430 

0.02711 

- 0.000 

11 

-0.06961 

-0.007 

0.01053 

0.01485 

- 0.000 

13 

-0.04357 

-0.004 

0.01048 

0.01265 

- 0.000 

15 

-0.02886 

-0.001 

0.01927 

0.01931 

0.000 

17 

-0.02001 

-0.001 

0.04803 

0.04780 

0.000 

19 

-0.01440 

-0.006 

0.16279 

0.16727 

-0.001 

21 

-0.01148 

-0.008 

0.64707 

0.65953 

-0.001 

Singularity 


CHEKX 

VAX 

VAX 

CHEKX 

Location 

dx 





1 

0.93197 

0.007 

-1.04498 

-1.05467 

0.001 

3 

0.93200 

0.007 

-1.04589 

-1.05567 

0.001 

5 

0.93208 

0.007 

-1.04870 

-1.05890 

0.001 

7 

0.93221 

0.009 

-1.05307 

-1.06466 

0.001 

9 

0.93240 

0.010 

-1.06097 

-1.07378 

0.001 

11 

0.93264 

0.006 

-1.07857 

-1.08830 

0.001 

13 

0.93293 

0.009 

-1.09704 

-1.10979 

0.001 

15 

0.93327 

0.006 

-1.13314 

-1.14448 

0.001 

17 

0.93366 

0.009 

-1.18730 

-1.20247 

0.001 

19 

. 0.93410 

0.009 

-1.29298 

-1.30950 


21 

0.93446 

0.006 

-1.38476 

-1 .40059 

0.001 
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Bcpl 
3r 

the normal velocity induced by all of the elements of the system 
(CHEKR or CHEKX) , the singularity-produced velocity (Vr or VAX), and the 
number of iterations required to achieve the no-flow condition. 

The singularity location numbering scheme used in the tables is 
as follows. The equally spaced nacelle singularities begin with point 1 
at the inlet plane and run axially to point 21 at the end cap plane. 
Similarly, the equally spaced end cap singularities begin with point 1 
on the nacelle centerline and run radially to point 21 at the nacelle 
surface. 

To establish the inlet flow field model, the nacelle generation 
procedure must reduce the normal velocities CHEKR and CHEKX along 
the nacelle and end cap, respectively, to an absolute value no greater 
than the no-flow criteria , EP. A comparison of these normal velocities 
to the selected value of EP (0.01) in Tables 1-5A clearly confirms the 
generality of this procedure with respect to the inlet-to-freestream 
velocity ratio (VELRAT) and the nacelle aspect ratio (AL'/2R'). 

The two cases presented in Tables 5A and 5B have identical 
dynamic and geometric conditions but differ by an order of magnitude in 
EP. Table 5B shows that increased accuracy is readily obtainable, but 
at the expense of additional iterations and consequently additional 


Bcpl\ 

Bx J ' 


velocity induced by the frees tream -inlet combination 


computer time . 
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At most of the singularity locations in Tables 1-5 the magnitude 
of the singularity-produced velocity (Vr or VAX) required to establish 

the no-flow condition differs substantially from that of the normal 

/ Sep 1 Sep 1 \ 

velocity induced by the frees tream-inlet combination or — J . 

This shows the importance of using a singularity strength adjustment 
scheme to generate accurately the inlet flow field model. 



aircraft speed, due to the larger size of the pre-entry streamtube. 

An additional pre-entry streamtube was calculated using a 5/1 
inlet-to-freestream velocity ratio but with a 2/1 nacelle aspect ratio. 
The points on this curve were indistinguishable from the 3/1 nacelle 



0 








© VELRAT * 60 
© VELRAT * 10 
(3) VELRAT = 5 
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aspect ratio case. This suggests that the path of the pre-entry 
streamtube is not strongly dependent on the nacelle shape. 

Table 6 presents the results of a check for continuity under- 
taken at three axial locations between several adjacent streamtubes 
in Figure 9 . At each of the locations , the discharge (V*A) was found 
by summing the V*A products of 20 subareas . The results of the check 
show the discharge to be nearly constant between streamtubes . The 
maximum variation of only 1 .2% from the average clearly demonstrates 
the precision of the streamline computational scheme. 


TABLE 6. RESULTS OF CONTINUITY CHECK 


AXIAL 

LOCATION 

DISCHARGE (V-A) 


Between 

Streamtubes 2 & 3 

Between 

Streamtubes 3 & 4 


-5.100 

2.57115 tt 

4.81250 tt 


3.000 

2.54358 TT 

4.79576 rr 


8.850 

2.51183 tt 

4.79334 tt 
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The Reingestion Example 

The dynamic and geometric conditions of the example problem 
are presented in Table 7. The example begins with the touchdown of a 
four-engined (wing-mounted) STOL transport and continues through the 
full deceleration process . 

The results are included in Table 7 and Figure 11. In this 
example, deceleration for the inboard engine occurs reingestion free. 
Also, the entrainment portions of both reversed jets never penetrate the 
pre-entry streamtubes of the engines discharging them. The aircraft 
configuration, however, proves to be highly prone to cross ingestion of 
the inboard engine exhaust to the outboard engine. Cross ingestion 
begins at an aircraft speed of about 70 miles per hour and continues 
through 50 miles per hour. Below this speed, the Maximum Penetration 
Point of the reversed jet of the inboard engine lies outside of the inlet 
flow field of the outboard engine. The entraining portion of the jet, 
however, continues to lie in this flow field and thus the possibility of 
further cross ingestion remains . 

Table 7 also lists the fluid particle time for those speeds 
where cross ingestion occurs. 
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TABLE 7. EXAMPLE PROBLEM 


DYNAMIC AND GEOMETRIC CONDITIONS 


VELJET 1 


— 

880 

Ft/Sec. 

INLET VELOCITY = 

440 

Ft/Sec. 

R 1 

= 

2.000 

Ft. 


AL' 

= 

22.6 

Ft. 


XJET' 

= 

20. 

Ft. 


YJET' 

= 

2. 

Ft. 


. ZJET' 

= 

0.0 

Ft. 


XSPACE' 

= 

16.0 

Ft. 


YSPACE' 

= 

24.0 

Ft. 


ZSPACE' 

= 

0.0 

Ft. 


DIAJET 1 

= 

2.30 

Ft. 


CIRCLE' 


0.0 

Ft. 


al 

= 

0.0° 



a2 


40.0° 




RESULTS 

Aircraft 

Runway 

Speed 

m .p .h . 

Inlet-to- 

Freestream 

Velocity 

Ratio 

Was 

Cross Ingestion 
Detected ? 

Fluid Particle 
Time (Sec.) 

90* 

3.33 

No 


80 

3.75 

No 


70 

4.30 

Yes 

.15903 

60 

5.00 

Yes 

.21740 

50 

6.00 

Yes 

.39996 

40 

7.50 

No 



^Touchdown speed 



FIGURE 11. THE EXAMPLE PROBLEM 




CHAPTER VII 


CONCLUSIONS 

This investigation succeeds in developing a method for analyzing 
the cross wind -free exhaust gas reingestion problem. The cases pre- 
sented pover a wide range of nacelle aspect ratios and inlet-to-freestream 
velocity ratios and clearly demonstrate the generality of the computer 
program . 

Results show the importance of using some type of singularity 
strength adjustment scheme in generating the inlet flow field model. At 
most points , the magnitude of the singularity-produced velocity required 
to establish the no-flow condition differs substantially from that of the 
normal velocity induced by the freestream-inlet combination. 

Additionally, data suggests that the shape of the pre-entry 
streamtube is uninfluenced by the nacelle aspect ratio. It appears that 
the accuracy of the method is independent of the nacelle shape, as is 
assumed in the development of the inlet flow field model. 
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APPENDIX I 


THE COMPUTER PROGRAM 


62 
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0001 

0002 

0003 

000 * 

0005 

000 * 

0007 

oooe 

0009 

0010 
oou 
0012 
0013 
001* 

0015 

0016 
0017 
001 * 

0019 

0020 
0021 
0022 

0023 

0024 

0025 

0026 
0027 
0026 

0029 

0030 

0031 

0032 

0033 

0034 

0035 

0036 

0037 
0036 

0039 

0040 

0041 

0042 

0043 

0044 

0045 

0046 

0047 

0048 

0049 

0050 

0051 
0092 


C 

C 

C 

C 

C 

C 


C 

c 

c 

c 

c 


c 

c 

c 

c 

c 

c 

c 

c 


c 

c 

c 

c 


c 

c 

c 


* STEP 1 * 


THE DIMENSIONING SECTION ***+*♦*+♦♦*♦** 

DIMENSION VR(Ji>#VAX<Jl}#0pR<Jl)/DpXCJl}4RADIAUtJl)*AXIAUCJl) 

DIMENSION VRI21 )/ VAXI21 )/ DPRI21 I/0PX121 1/RADl AL (21 )/ AXI AL (21 ) 
DIMENSION K0H4 1 21 1/KUN3 ( 21 ) 

DIMENSION CHEKR(21)/CHEKX(2U 


THE INPUT SECTION 
PROGRAM UPTIUNS 

Hast, nojeT/nocros/Nspeeo# nocard/ nopest) 


LAST«2 

N0JET*0 

NOCR0S»l 

NSPEE0*0 

NOCAKO'O 

NOPeSTpO 


PROORAM INPUTS 
(EPiDEkTAX) 

E |NPVT DEITAX AS POSITIVE FOR FORWARD PLOTTING/ NEGATIVE FDR 
REVERSE. 

DELTAX«Ot 300000000000000000000 


DYNAMIC INPUT 

(U8A/ VELRAT/ VELJET) 

UBA-88, 

VELRAT«5, 

VELJET* 8 80.00000000000000000000000000 


GE0METR1 | R jJj^J X j eT# YjgT/2jET/ALPHAi/ALPHA2/XSPACE/YSPACE/ISPACE/ 
DJAJET/CIRCLE ) 

M2. 000000000000000 

ALA*22. 600000000000000000000000 

XJET* 20, 000000000000000000000000 

Y JET* 2. 00000000000000000000000000 

IJET«0. 0000000 

ALPHA1 *7. 50000000000000000 

ALPHA 2 *40.000000000000000000000000000 

XSPACE* 16, 000000000000000000 

YSP ACE *24 ,000000000000000000 

2 SPACE *0,00000000000000000000000000 

D! A JET*2, 30000000000000000000000 

C I RCLE'0, 0000000000000000000 


PROORAM ADJUSTMENTS 

IXCHEX/ SLOPE/ J6/J1/KLAST/D1V2/DIV4/N3ST0P) 


I 
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0993 

0094 
0099 
0096 
009? 

0098 

0099 
0060 
0061 
0062 
0063 
0066 
0069 
0066 
0067 
0069 

0069 

0070 

0071 

0072 

0073 
0076 
0079 

0076 

0077 

0078 

0079 

0080 
0081 
0082 
0083 
0086 

0085 

0086 

0087 

0088 

0089 

0090 

0091 

0092 

0093 
0096 

0095 

0096 

0097 

0098 

0099 

0100 
0101 
0102 
0103 
0106 

0105 

0106 

0107 

0108 

0109 

0110 
0111 
0112 
0113 
0116 

0115 

0116 
0117 

,0118 


XCHEK.50, 

SLOPEi. 75 
J6.5 
Jl*21 
KLAST«20 

DIV2.6.000000000000000000000 
01 V6"2. 

Nt$TQP"6 

C *****♦*•*♦• 

c * STEP 2 * 

C »••*****«** 

C • INITIALIZATION I 
OHOLOvOELTAX 
0Sl-2,*VELftAT*U8A 
NSP.O 

VELCET«V£LJET 
JL AST* J 1-1 
01 VI" JLAST 
AL"Al A/R 
0X-AL/01V1 
DSR-i./DJVl 
UAST.JLAST 
0T1.6.28318/01V1 
I1-1LAST*! 

X7-ILAST 
K8"tl 
NORAD. 1 
N0AX«1 
OFF.O, 

KTEST"6*J1 
7011 CONTINUE 

0P5.l./<4,*3,16199*u8A) 

K0N1*0 
K0N5«1 
OUTPUT K0N5 
C 

IF (NOCARD. EQ.O )60 TO 2 
RE AOI 5/ 396 )VR 
READ! 5/ 396) VAX 
00 7008 J.l/Jl 

7008 WR1TE{6,383)VR< J)*VAX(J> 

CO TO 7010 

C ••*******♦• 

c • STEP 3 • 

C *•♦*♦****♦* 

2 KON 1 "KON 1* 1 

IF(KONi,0T.KLAST)6D TO 1007 

!F(LA5T.EQ.1)WRITE(6«306)K0N1 

JF(LAST,E0,l)WRlTE(6/300) 

IP(LAST.EQ.1)WR!TE<6*305) 

IF ( LAST. EQ. 1 (WRITE (6*300) 

C 

C INITIAL 11ATI0N II +*«+♦+*♦+*«♦++♦♦*+++++*+++++♦+++«**++++++ 

K0NX»0 
KONR.O 
DO 3 J"1/J1 
AJ-J-l 
X.AJ*DX 
SR.AJ6DSR 

IF( J,CE,Ji)X«X-0X*.5 
1F(J,0E.JV)SR«SR-0SR*,5 
IF ( KUN1 .CT .1)00 TO 6 
C 

CALL COOP(X*5R/KQN5,A6*B6/A3/B3/SUH2/SUM3*$UM6>SUM5*0Tl*0P16/Tl. 

t RAD.N0AX.il/AL/N2.N6) 

AXJAL(J)"l,-QSl*0P5*AL'*SUM3 


\ 



oil* 




RADIAL ( J ) »-OSl*OP5*SOM2 

0120 




VAX t J ) AXIAU J) 

0121 




VR( J I ••RAD J At, ( J ) 

0122 




IFUXIAL(J) 15,29/29 

012) 


5 


*6176(6,307) 

0124 




00 TO 1007 

0125 


29 


IF (LAST, EO. 1 ) WRITE (6,306 ) J, OFF, RADIAL N), OFF, OFF, OPP/ AXIAL (J>< 01 

0126 



t 

off,qff,uff 

0127 




00 TO 3 

0126 

C 




0129 

c 



* STEP 4 * 

01)0 

c 



*********** 

01)1 


4 


CALL JTUFF<J,X,SR,VAX,VA/AL,0X,DSR,DT1,K7>K6,KON5,TRH33,TRM34,TI 

01)2 



t 

5,TRM36,NORAD,NOAX,NZ,N4) 

01)) 




TRMJ6.0, 

01)4 




0PR<J)«TRM34*TRM35 

01)5 




DPXt J)«TRH33*TRH36 

01)6 




CHEKR(J)»DPR(J)*VR( J)*RADIAL(J> 

01)7 




CHEKX(J)»OPX(J)*VAX(J)+AXIAL(JJ 

01)6 

c 




01)9 

c 



* STEP 5 • 

0140 

c 




0141 




IF(A65(CHEKX(jn.LE.EP)GO TO 43 

0142 


S3 


KONA ( J ) ■ 1 

014) 




GO TO 44 

0144 


43 


KONA ( J ) *2 

0145 


44 


IF(A6S(CHEKR(jn.LE.EP)C0 TO 4) 

0146 


)1 


K0N3 ( J ) *1) 00 TO 46 

0147 


45 


K0N3 ( J ) »2 

0146 


46 


KONR»KQNR*KON)(J) 

0149 




K0NX«K0NX*K0N4(J) 

01)0 




IF (l AST , EQ , 2 ) 00 TO 3 

01)1 




KR1TE(6,306)J, VR( J), RADIAL ( J J* TRHB4,TRH35,VAX( J >, AXIAL! J >,TRM33, 

01)2 



t 

K36,CHEKRN),CHEKX(J) 

0153 


1 


CONTINUE 

0154 




KSUM*KONX«KONR 

0155 




I F ( LAST, EQi 2 )G0 TO 6000 

01)6 




*RITE(6,300) 

0157 




WRIT6(6,301)KSUH 

01)6 




WRITE (6,303 )KONR 

01)9 




WRITE (6, 302 IKONX 

0160 




WRITE ( 6, 300 ) 

0161 




WRITE (6, 300 ) 

0162 

8000 


CONTINUE 

016) 




IP(K5UM,CE,KTEST)C0 TO 26 

0164 




IF ( K0N1 a EQa 1 ) 60 TO 2 

0165 

C 




0166 

C 



• STEP 6 * 

0167 

c 




0166 




CALL STABLE (Jl, KONA, K0N3,eHEKX,CHEKR, VAX, VR) 

0169 




00 TO 2 

0170 

c 




0171 


26 


WRITE (6, 313 ) K0N1 

0172 




IF ( LAST a OE ■ 2 } 00 TO 7010 

0173 




CO TO 7007 

0174 

c 




0175 

c 



• STEP 7 * 

0176 

c 




0177 

7010 


CONTINUE 

0178 




IP ( NOPEST aOT » 0 ) 00 TO 1602 

0179 




IP(NUJET.EQ.O,AND.NOCROS,EQ.O)GO TO 1006 

0160 

1602 


CONTINUE 

0161 




VELJET.VELOET 

0162 




CALL THE JET (NOCROS/YJET,ZJET, ALPHA 1,ALPHA2,DIAJET,VELJET,XJET,XM 

0163 



t 

,RMPP,2SPACE,YSPACE,XSPACE,XOP,RCROSS,U8A,R,THEATA,NSP,RJET) 

0164 

C 






0119 

1006 

KQN9«2 

0166 


OUTPUT K0N9 

0167 


N1»0 

0168 


K00P1*1 

0169 


Kft ITE ( 6/ 314 ) 

0190 


NOP£»i 

0191 


NAL-0 

0192 


TIM£1«0. 

0192 


IP (NOJET « EQ , l )G0 TO 1437 

0196 


IP (NQCROS |GT i 0 ) GO TO 1600 

0199 


1P(NUPES7,0T,0)C0 TO 1437 

0196 

1603 

CONTINUE 

0197 


READ15<306)X/SRjTIME1 

0196 


00 TO 1436 

0199 

1417 

X»XMPP 

0200 


SR. RUPP 

0201 


WRITE (6/ 397 )Xi SR 

0202 


SR-SR-CIRCLE/2, 

0202 


SlOP62«<RJ£7-<OUJET/2.>*COStTMEATA)-SR)/(XJET-X) 

0204 


BEE.4R-SL0PE24X 

0209 


IP<NOPEST.0T,O)6O TO 1603 

0206 


00 TU 1601 

0207 

1600 

X«XOP 

0208 


SRaRCROSS 

0209 

1601 

CONTINUE 

0210 


NAL«A8S { X ) / AL 

0211 


SPINAL. OT. 0)00 TO 1439 

0212 


NAL»1 

0212 

1419 

RNAL'NAL 

0214 


0ELTAX.(AB5(X) )/(OIVi*RNAU 

0219 

C 


0216 

C 


0217 

c 


0216 

1436 

CONTINUE 

0219 


WRlTE(6/308) X/SRyXCHEK 

0220 


K00P2»XCHEK 

0221 


N2»l 

0222 


N3«0 

0222 


N4«l 

0224 


LESSul 

0229 


XXi. 0,0000000000 

0226 


XX2.0X 

0227 


NOOO* 1 

0226 


DDX-0.000000000000000000000 

0229 


J9-0 

0220 


NTIM»0 

0231 


n:«i 

0232 


NSSUH.O 

0233 


NZEA«1 

0234 


EP2»AB5 ( DELTAX ) 

0239 


EP3— EP2/3, 

0236 


EP4«Al*EP2«EP3 

0237 


J3MOOOO 

0236 


IP(DELTAX)1022/1023/1023 

0239 

1022 

J»J 1*2 

0240 


J2.-2 

0241 


l P ( X . LT , EP4 ) LESS* 2 

0242 


00 TO 1006 

0243 

1023 

J»-l 1 J2»2 

0244 


JP(X,GT,EP3)NOPE«2 

0249 

1006 

IP(OKLTAX) 1040/ 1041/1041 

0246 

1040 

JP(J3,E0,J)G0 TO 1020 

0247 


1P(X.LE,EP4, AND, LESS. EQ. 1)00 TO 1042 

0246 


00 TU 1026 

0249 

1042 

J3.JJG0 TO 1028 

0290 

1041 

1 P ( X. OE , EP3 ) 00 TO 1020 


• STEP B * 



67 


0291 



00 TO 1028 

0292 

C 



0299 

1020 


N2-2 

029* 



IF(J2>1029, 1026, 1026 

0299 

1029 


X«AL-,5*0XJ JST0P-1JDELTAX«-2,*DX 

0296 



J 3- 1 000 

0297 



00 TO 102* 

0298 

C 



0299 

1026 


IF ( NOPE ,EQ , 2 ) 00 TO 1902 

0260 



X-0.0000000000000000000000000 

0261 

1902 


JSTOP- J1 

0262 



J3-J1-2 

0263 



J«OjNO:V2iD!V2jND!V3«01V2/2. 

026* 



0ELTAX»2t*DX/DIV2 

0269 



NQGO-2 

0266 



J2B-1 

0267 



00 TO 1069 

0268 

C 



0269 

102* 


J-J-J2 

0270 



IF ( J 1 80 1 JSTQF ICO TO 10*9 

0271 



CO TO 1028 

0272 

C 



0273 

10*9 


N2-3 

027* 



I P ( J i EQ , 1 ) GO TO 1027 

0279 



X«AL-,5*DX 

0276 



J3-1000 

0277 



RN35TP-N35T0P-2 

0278 



OELTAX-OHOlD/RNSSTP 

0279 



N2-* 

0280 



00 TO 1028 

0281 

C 



0282 

1027 


DELTAX-OHOLO 

0283 



J«0 

028* 



X«0. 

0289 

C 



0286 

1028 


IPISR.LE.O. ) WR I TE { 6, 386 ) 

0287 



IF< J,EQ.l, AND. SR. LE, 1,00)00 TO 1203 

0288 



CALI C00P<X,SR,KUN5,A*,BA,A3,B3,SUM2,SUM3,SUM*,SUM3#DTl,0Pi*,Tlj 

0289 


c 

RAD/NOAX/ 1 1, AL/NZ/N*) 

0290 



I F ( NOR AD i E Q ( 2 ) GO TO 10*8 

0291 



DPOTR— Q51*OP9*SUM* 

0292 

10*8 


1 F (NOAXa EO , 2 ) GO TO 1021 

0293 



OPOTX*l.-QSi*OP3*X*SUM3 

029* 

C 



0299 

1021 


CALI 5TUFF! J,X,SR, VAX,VR, AL,DX,DSR,DTl,K7,KB,K0N5,TRH33,TRM3*,Tf 

0296 


t 

5,TRM36,NORAO,NOAX/NZ,N*) 

0297 



IF (NCJRAO ■ EQ • 2 ) GO TO 10*9 

0298 



VRAD-0P0TR+TRM3A+TRM39 

0299 



1F(NUGO.GT.2,ANO.VRAO.GE.O. 1VRAD.0. 

0300 

10*9 


IF (NUAX i EQ • 2 ) GO TO 1090 

0301 



VAXY«0P0TX+TRM33*TRM36 

0302 

1030 


GO TO < 1031, 1091, 1066, 1067, 1076 ),N* 

0303 

1076 


NORAD-1 

030* 

1051 


IF(VAXY,l£.0))WRJTE(6,387) 

0309 



OROX-VRAO/VAXY 

0306 



IF(ABS(VRAO)«LE,EP ) DRDX-O, 

0307 



I F ( NOPE f EQ 1 1 ) GO TO 1500 

0308 

C 



0309 



1 F ( VAXY , GE aO ' I GO TO 1301 

0310 



NOPE -3 

0311 



00 TO 1101 

0312 

1901 


NOPE-i 

0313 

1900 


CONTINUE 

031* 

C 



0319 

C 


0316 



I F ( NUGQ , t Q i 1 ) GO TO 1096 
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i 


0117 

0 * 1 # 

0*19 

0*20 

0*21 

0*22 

0 * 2 * 

0*2* 

0*29 

0*26 

0*27 

0 * 2 # 

0*29 

0*30 

0*31 

0*32 

0*3* 

0334 

0339 

0336 

0337 
0*3# 
0*39 

0340 

0341 

0342 
0*43 
0344 
0349 
0*46 
0*47 
034# 
0*49 
0*90 
0*91 
0*92 
0*93 
0394 
0399 
0*96 
0*97 
0*9# 
0399 
0360 
0*61 
0362 
0*63 
0364 
0369 
0*66 
0*67 
0*68 
0369 
0*70 
0*71 
0*72 
0*7* 
0*74 
037* 
0*76 
0*77 
0*7# 
0*79 
0*60 

0381 

0382 


14)0 


1402 

1403 

1404 
1409 

1406 

1407 
14*2 


1429 

1433 

14)6 


1431 


14*4 

14)3 


1096 


l)(ABS(0RDX),LTiSL0PE)NSL0PE«l 

IF(A8S(OROX),OE. SLOPE INSL0Pe«2 

I F ( NS SUM , E Q , 0 ) GO TO 14*0 

IFINZER. £0,2. AND. Nl.GT, 0)60 to 1096 

lPtNl,EQ,0)SU TO 1432 

N)SUM«NSSUM+NSLOPE 

00 TO 1096 

NSHOIO'NSLQPE 

NSSUM.NS10PE 

JHQlD* J 

XX1H«XX1 

XX2H»XX2 

VRHOIO«VRAO 

VXH010>VAXY 

OROXHmORDX 

NUH#ER«1 

SO TO (1096* 1402/1403; 14 04* 1403* 1406/ 1407)*N0GQ 
lF(ABS(DROX)|GT a SLOPE)GO TO 1408 
IF( Jt G£<tl6)GQ TO 1077 
00 TO 1096 

1F(AS3(DR0X) « LT, SLOPE ) GO TO 1082 
00 TU 1096 

if ( abs < dRox ) «gt. slope too to 140# 

00 TO 1096 

IF(ABS(OROX).LT.SlOFE)CO TO 1200 
GO TO 1096 

]F(ASS(OROX) l GT.St.OPE)GO TO 1409 
00 TO 1096 

1 F ( ABS ( OROX) (GT (SLOPE ) GO TO 1*09 
CO TO 1096 

IF(N6H0L0<EQ«2)GO TO 1429 

1F(N3)UH,LE.6)G0 TO 14)0 

X'XHOLO 

SR*SRHQ10 

NSHOIO'2 

CO TO 1431 

IF(N5SUM,G£,6)00 TO 14*0 

00 TU 1436 

NUHBER»0 

NZER»1 

X»SRM010 

5R.XH0L0 

NSHOtO* 1 

VRAO'VRHQIO 

VAXY»VXHOLO 

DADX«DR0XH 

N1 »0 

NSSUM«NSH010 
NUMBER«NUMBER*1 
I F (NUMBER I GE , * ) GO TO 1433 
WRITE (6* 399 ) 

WRITE (6* 300 ) 

J P ( J i E9 1 JHOLD) GO TO 1434 

XX1-XX1H 

XX2»XX2H 

J» JHOLD 

00 TU( 1096* 1408/ 1082* 1408* 1200* 1409* 1409) *NOCO 
WRITE (6* 398 ) 

00 TU 1007 

CONTINUE 

call RUNOE ( X* SR* XMOLO/SRHOIO/OROX/OEITAX* J3* At* EK1* EK2* EK3*EK4,N 
J/OX/OOX* J3*N0G0*XXl*VAXY, VRAO, EP/NZER ) 

I F INZER • EG i 2 ) GO TO 1439 

00 TO (10 1 2* 107 5* 1075* 1013 )*NJ 

00 TU ( 1028* 1073* 1066* 1080* 1103* 1080*1073 )*N000 


1073 



om 


00 TO 1028 

0114 

1012 

CONTINUE 

0)19 


GO T0( 1020/ 1073/ 1006/1080/1103/1080/ 1073I/N0C0 

0306 


00 TO 1020 

0307 

101) 

WRITE (6/ 383) X/ SR 

0300 


WRZTE(6/900) 

0309 


Ni»0 

0)90 

C 


0391 


IP(NQJ£T,GT.O,OR.NOCROS,GT.O)6Q TO 1613 

0392 


10 (NUftEST.EQ i 0 ) GU TO 1097 

0)9) 


1P<NAL,EQ.-10)GO TO 7009 

0)96 


BORDEA«SL0PE2»X*8EE 

0)99 


1P(SR,GT,BURD£R)NAL—10 

0396 


IPINAL.£Q.-10)WRITE(6/393) 

0397 

1613 

CALL TIMEiNT!H/XA/XB/RA/RB/VRA/VPB/X/Sft/VAXV/VRA0/T!HEl/T!HE2) 

0)96 

7009 

CONTINUE 

0399 

C 


0600 

1097 

00 T0( 1003/ 1062/ 1006/ 1062/ 1103/ 1062/ 1062 )/N000 

0601 

C 


0602 

1003 

K00P1«K00P1«1 

0603 


IP (K00P1-K00P2 11029/ 1029/7007 

0606 

1029 

GO TU( 1006/ 1026/ 1028/ 1066 )/N2 

0609 

1066 

N3.N3*1 

0606 


IP(N3,LT,N3STOP)GO TO 1020 

0607 


oeltax«ohqlo 

0600 


N2«3 

0609 


00 TO 1020 

0610 

C 


0611 

1069 

J6-0 

0612 

1061 

J5«0 

0613 


J«J*J2B 

0616 


1PU20.lt. 1)G0 TO 1100 

0613 


IPU.EQ.1IG0 TO 1096 

0616 


XX1«XX2 

0617 


IPU.EO.l ) XX 1*0. 0000000000 

0610 


XX2»XX1*DX 

0619 

1096 

I P ( NOGO , EO, 3 )G0 TO 1086 

0620 


00 TU 1076 

0621 

1100 

XX2/XXI 

0622 


XX1.XX2-0X 

0623 


IPU.EO, 2 )XX1»0, 000000000000000000 

0626 


IPUiLT.l. AND. Nl.EQ. 0)G0 TO 1203 

0629 


1 ft (NQGQ, E9, 3 )G0 TO 1086 

0626 

1076 

N6«2 

0627 


I F (NOPE . EQ. 1 ) GO TO 1028 

062S 


NTRY»2 

0629 


00 TO 1087 

0630 

1062 

J6« J6*l 

0631 


IPU6.E0.N0IV3. AND. NOGO. E0.4)G0 TO 1081 

0632 


JPU6,EQ,N0IV2. AN0.NUGU.EQ.6IG0 TO 1081 

0633 

1202 

IPU6.E0.N0IV2, AND. NOGO. EQ.61G0 TO 1201 

0636 


IPU6.EQ.N0IV3.AND. NOGO . E 0.6) GO TO 1201 

0635 

1092 

IPU4.E«iN0IV3)GU TO 1061 

0636 

. 

1 p U6.EQ.N0I V2 )G0 TO 1065 

06)7 

1073 

IP U3.EQ.ND1 V2 )G0 TO 1076 

0638 


GO TO 1080 

06)9 

C 


0660 

1077 

N0G0*1 

0661 


N6*l 

0662 


DELTAX«2«*0X 

0663 


00 TO 1096 

0666 

C 


0669 

1001 

N0G0*2 

0666 


0ELTAX»2,*DX/0IV2 

0667 


J2B«1 

0668 


CO TO 1092 



04*9 

0490 

0491 

0492 
0499 
0494 
0499 

0496 

0497 

0496 
0499 

0460 

0461 

0462 
0469 
0464 
0469 
0466 
0*67 
0*68 

0469 

0470 
0*71 
0472 
0479 
0*74 
0479 
0*76 
0*77 

0478 

0479 

0480 

0481 
0462 
0489 
0*84 
0489 
0486 
0*87 

0488 

0489 
0*90 
0*91 
0492 
0*99 
0494 
0*99 
0*96 

0497 

0498 
0*99 

0900 

0901 

0902 

0903 
090* 
0909 

0906 

0907 

0908 

0909 

0910 

0911 

0912 

0913 

0914 


C 

1408 NTRY*1 

1087 N0G0*3 

XMAX».999*X 

J2B.1 

I F ( J4 1 IE t NDI VI ) CD TO 1093 
1109 J4«N0!V3 

061200—1. 

00 TO 1106 
1093 J4*0 

081200*1. 

C NeXT CARO gets a DELTA R FOR THE FLIPPED runge-kutta method 

1106 DELTAX—0X/(D1V4«D1V2) 

] F (NTRY.EQ. 1 )GQ TO 1410 
:F(NOFE.E0.9)GO TO 1901 
] F (NUGO. EO. 9 ) GO TO 1103 
GO TO 1086 

1084 J4.J4*0F1200*ND1V3 

0P1200—0P1200 
| F 1 NOGO , E Q , 9 ) CD TO 1103 
1086 J F ( X.OE ,XX2 )G0 TO 1061 

IF(X«1T.XMAX,ANO,N1.EQ.O)GO TO 1428 
1429 XMAX*X 

GO TO 1107 

1427 JF(VAXY)1426<1426#1087 

1428 IF (VAX Y >1101/142 9.1429 

1109 IF(X.LE,XX1)G0 TO 1061 

IF (X. GT ,XMIN . AN0.N1 , EQ,0 100 TO 1427 
1426 XMJN«X 

1107 IF (SR ,GE . 1 . 01 )G0 TO 1080 
WRIT£(6,3B8) 

GO TO 7007 
C 

1082 NOG0*4 

DELTAX«2,*(XX2-XJ/0IV2 

J28*i 

J9«0 

GO TO 1096 
C 

1409 NTRY'l 

1101 NOGO-9 

XMIN*l i001*X 
J28.-I 

JF(J4|GE.NDIV3)G0 TO 1093 
GO TO 1109 
C 

1200 N0G0"6 
DELTAX»«2,*(X»XX1 } /DIV2 
J9*0 

J28— 1 
GO TO 1096 
C 

1201 N0G0*7 
DEITAX*-2.*0X/0IV2 
J2B.-1 

GO TO 1092 
C 

1080 N4*3 

N0AX*2 

XX3"X 

X*XX1 

GO TO 1028 
1066 N4*4 

X*XX2 
VRi-VRAD 

GO TO 1028 . . . ... 
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0519 

1067 

0516 


0517 


0518 


0519 


0520 


0521 


0522 

C 

0523 

1*10 

0526 


0525 

C 

0526 

C 

0527 

C 

0528 

1203 

0529 


0530 


0931 


0532 


0533 


053* 

C 

0535 

C 

0536 

c 

0537 

7007 

0938 


0939 


0540 


05*1 


0962 

1007 

05*3 

300 

05** 

301 

05*5 

302 

09*6 

303 

05*7 

30* 

05*8 

305 

05*9 


0550 


0951 

306 

0552 

307 

0553 

308 

055* 

313 

0555 


0556 

31* 

0557 

383 

0558 

386 

0559 

387 

0560 

388 

0561 

390 

0562 

391 

0563 


056* 

392 

0565 


0566 

39* 

0567 

395 

0568 


0569 

396 

0570 

397 

0571 


0572 

398 

0573 

399 

057* 


0579 



N*«5 

X-XX3 

VR2-VRAD 

VRAD" ( ( DX-DDX) /OX )* ( VA1-VR2 >*VR2 

NOAX«l 

NORAD-2 

CO TO 1028 

NTRY-2 
00 Tg 1096 

• STEP 9 * 


WRITE ( 6/ 390) 

TIME1-T1ME1*R/U8A 
IF(NOJET,GT,0)WRITE(6/391)TIME1 
IF (NUCROS, GT.O) WRITE (6/392 JTIHE1 
I F (NO JET, EQ.O.AND.NOCROS.EQ.O) WRITE (6/39*1 
CO TU 7007 

• STEP 10 * 


VP RADIAL VR-RAD 

VAX-RAD VAX-AX CHEKA 


V 

CHE 


{FtNSPEEO.EO.OIGO TO 1007 
NSP-NSP+1 

1P(NSP,EQ.NSPEE0)G0 TO 1007 
R8A0!5/396)U8A 
CO TO 7011 
STOP 

FORMAT! 'O' 1 
FORMAT! 'KSUM»</ 15) 

FORMAT ! ' KONX- 1/15) 

F0RMATCK0NR-I/I5) 

FORMAT < ' K0N1- '/IS) 

FORMAT! ' J 

l AX VAX AXIAl 

Cl) 

F0RMAT(2X/I10/8(2X/F10,5)/2{2X/F6,3>) 

FORMAT! 'THE NACELLE CAN NOT BE GENERATED!) 

FORMATOFIO.5) 

FORMAT! 'THE INLET FlOW FIELD MODEL HAS BEEN GENERATED- THE NUMB 
t OF ITERATIONS REQUIRED WAS'/I5> 

FORMAT! ' X SR XCHEK') 

FORMAT ( 8F10, 5 ) 

FORMAT! 'THE RADIUS HAS BECOME NEGATIVE') 

FORMAT! 'THE AXIAL VELOCITY HERE 15 NEGATIVE') 

FORMAT ( ' THE RADIUS HAS PENETRATEO THE NACELLE') 

FORMAT! 'THE STREAMLINE HAS ENTERED THE INLET') 

FORMAT! 'EXHAUST GAS RE-IN JESTED '/ FB. 5/ 'SECONDS AFTER PENETP 
C ION OF THE PRE-ENTRY STREAM TUBED 

FORMAT! 'EXHAUST GAS CROSS- INJESTEO ' , FB . 5/ 'SECONDS AFTER PEN 

t RATIUN OF THE PRE-ENTRY STREAM TUBED 

FORMAT! 'EXHAUST GAS INJESTION WA5 NOT DETECTED') 

FORMAT! 'SECTIONS OF THE ENTRAINMENT PORTION OP THE JET LIE WITH! 
t THE PRE-ENTRY STREAM TUBED 
FORMAT ( 8F 10 ■ 5 ) 

FORMAT! 'THE MAXIMUM PENETRATION POINT OCCURS AT'/F10.5/ ' 

C D 1 1 AXIALLY AND'/F10,5/ 'RADII RADIALLY' > 

FORMAT! 'THE POINT IS CURRENTLY UNOBTAINABLE'! 

FORMAT! 'PROCEDURE CHANGED AT THIS POINT, ABOVE NUMBERS ARE NO C 
t O') 

END 
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0001 

0002 

0003 

000 * 

0005 

0006 

0007 

0008 

0009 

0010 
0011 
0012 
0013 
0016 

0015 

0016 
0017 
0016 

0019 

0020 
0021 
0022 


0001 

0002 

0003 

0006 

0005 

0006 
0007 
0006 

0009 

0010 
0011 
0012 
0013 
0016 

0015 

0016 

0017 

0018 

0019 

0020 
0021 
0022 
0023 
0026 

0025 

0026 

0027 

0028 

0029 

0030 

0031 

0032 

0033 
0036 

0035 

0036 

0037 

0038 

0039 


SUBROUTINE T|ME(NTIM/XA/XB/RA/RB/VPA/VPB/X/SA/VAXY/VRAD/T1ME1/ TI 
8 2 ) 

C THIS subroutine calculares dimensionless time 
JP(NT1M, £0.0)00 TO 1 
XA«XB 
RA»R8 
VPA*VPB 

1 XB.X 

RB.SR ! 

VPB«5QRT(VAXY**2*VRA0**2) 

| P (NTIM, EQ| 0 ) 00 TO 2 
00 TO 3 

2 NT JM«1 

00 TO 5000 

} VAB.(VPA*VPB)/2. 

XAB-XA-XB 

RAB-RA-RB 

DSAB*SQRT ( XAB**2*RAB**2) 

TAB»DSAB/VAB 
T!ME1»T1ME1*TAB 
5000 RETURN 

ENO 


SUBROUTINE RUNGS <X/SR/XHDLD,SRHDID,DRDX,DELT AX/ J3, AL/EK1/EK2/EK- 
8 K6/N1/J/0X/0DX/ J3/NOGO/XX1/VAXY/VKAO/EP/N2EA) 

C THIS SECTION CONTAINS THE RUNGE-KUTTA METHOD, 

Nl*Ni*l 

IP(NUGO,EQ,3,OR.NOSO,EQ.5)GO TO 6 
00 TO 1079 

6 N5«l 

IP(A8S!VRAO),LE,EP)GO TO 7 

drox»vaxy/vrao 

IP(ABS(VAXY)|LE.EP)0RDX>0. 

00 TO 2 

7 NZER»2 

GO TO 5000 

5 N5»3 

GO TO 2 

6 J5» J3+2 

3 DDXaX-XXl 

GO TO 5000 

1078 N5»2 

2 XSWAPaX 

RSNAPaSR 
XaRSNAP 
SR.XSWAP 

GO TU( 1079/ 3000/31/N5 

1079 GO T0( 1016/ 1017/ 1018/ 1019I/N1 

1016 EK1.0RDX*DELTAX 
IP(J3,E0, 1000)00 TO 1063 

1066 XHOLD-X 

SRH0L0»5R 

X»X*0ElTAX/2, 

SR»SR*EKl/2, 

GO T0(5000/3/5/3/5/3/3)/N0GO 
1063 X*ALI J3*10000j GO TO 1066 

1017 EK2*0R0X*0ELTAX 
SR»SRH0lD*EK2/2. 

GO 70(5000/5000/ 1078/5000/ 1078/ 3000/ 5000 )/NOGO 

1018 EX3>DR0X*DELTAX 
X«XHOLO*OELTAX 

IP(J.E0,J3)X«AL-,5*OX . 
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0040 

0041 

0042 

0043 

0044 
0043 
004* 
0047 
004* 


SR*SRH0L0*£K3 

00 TQ<3000,4,3,4/5/4/4)/N0G0 
1019 EK4*OAOX*OELTAX 

DElRA0«(EKl*2,*EK2+2,*EK3*£K4)/6, 

X-XHQLD+DELTAX 

5R*SRHOL0*DELRAD 

00 TO <3000/ 3000/ 1078/ 3000,107*/ 3000/ 3000 )/NQGO 
3000 RETURN 

END 


0001 

0002 

0003 

0004 
0003 
0006 
0007 
000 * 

0009 

0010 
0011 
0012 
0012 

0014 

0015 

0016 
0017 
001 * 

0019 

0020 
0021 
0022 

0023 

0024 

0025 

0026 
0027 
002 * 

0029 

0030 

0031 

0032 

0033 

0034 
0033 

0036 

0037 

0038 
00*9 
0040 


SUBROUTINE THBJET<N0CR0S,YJET,ZJET/ALPHAI, ALPHA2,DIAJET/VELJET/X 
t T/XMPP/RMPP/ZSPACE/ YS PACE/ XS PACE/ X0P/RCR0SS/U8A/R/ THE ATA/NSP/RJE 
C THIS SECTION CONTAINS THE LOCKHEEO CORRELATION, 

IP {NSP i CT , 0 ) 00 TO 200 
XJET«XJET/A 
YJET»YJET/R 
ZJET»ZJET/A 
0JAJET»DI AJET/R 
At PHAl»ALPHAl/37 ,29378 
ALPHA2«ALPHA2/ 37, 29378 
200 CONTINUE 

VELJET«VE1JET/U8A 

AJET»SQAT<YJET**2*ZJET**2) 

IPIALPHAZ.LT. ,01 ICO TO 1 

BETA* ATANI TAN <ALPHA1)/S!NIALPHA2>> 

IPIALPHAl.LT, ,01)8ETA*0,000000000000000 
GO TO 2 

1 BETA*90. 000/37, 2937* 

2 continue 

TH£6TAkAC0S<C0S<ALPHAU*C0S<ALPHAZ>) 

A*<1,.,734*<SJN<THEATA) >**.683) 

PHPP«2, 97*01 AJET*A*(VEL JET**, 94> 

qmpp*pmpp*tan<theata> 

XHPP»XJET-PMPP 

RMPP*SQRT(RJET**2*QHPP**2*2,*QMPP*(YJET*C0S(BETA)*ZJ6T*SIN(BETA) 
OUTPUT XHPP/RMPP 
1P(NOCROS,EQ.O)GO TO 3000 
I P ( NSP , GT . 0 ) GO TO 203 
XSPACE-XSPACE/R 
YSPACE-YSPACE/R 
ZSPACE»ZSPACE/R 
203 CONTINUE 

ZMPP*RMPP*SIN(8ETA> 

ZOP«ZSPACE*ZHPP 

YMPP*RMPP*COS(BETA> 

YOP.YSPACE-YMPP 
XOP«XMPP-XSPACE 
RCR0SS«SQRT(Y0f , ****Z°9**ZI 
3000 RETURN 

END 


0001 

0002 

0003 

000 * 

0003 

0006 

0007 


SUBROUTINE STABLE I Jl, KONA, K0N3/CHEKX/ CHEKA, VAX/ VR ) 

DIMENSION K0N4< J l ) , K0N3 < Jl ),CH£KX< J1 >, CHEKA < J l >, VAX < J 1), VR < J 1 ) 
C THIS SECTION ADJUSTS THE COMPENSATORY SINGULARITY STRENGTHS, 

40 00 *2 I « 1/ J 1 

1PIK0NA! J).GE,2)G0 TO 47 
A102*CHEKX( 11/2. 

1P(VAXU)*CHEKX(U 164/64/62 
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0008 

62 

0009 

63 

0010 

66 

0011 

0012 

67 

0013 


0016 


0013 

86 

0016 

60 

0017 


0018 

61 

0019 

62 

0020 


0021 



IF <A8SiVAX(I},-ABS(Al02» ,63/63/66 
VAX( | J ■ VAX ( I) /2 >09 
CO TO 67 

VAX<i>-VAX<l)-A102 

IF,KUN3< J ),GE,2)C0 TO 62 

A101»CH£KR,I)/2a 

IP(VR, I , 6CHEKR, I) >61/61/56 

JHABSIVRI I> 1-ABStAlOU >60/60/61 

VA<,»*VRf,)/2.0S 

CO TO 62 

VR, 1 , «VR, I > — A 1 0 1 

CONTINUE 

RETURN 

END 


0001 

0002 

0003 

000 * 

0009 

0006 

0007 

0006 

0009 

0010 
0011 
0012 
0013 
0016 
0013 
0016 
0017 
00i« 

0019 

0020 
0021 
0022 
0026 

0025 

0026 
0027 
0026 

0029 

0030 

0031 

0032 

0033 
0036 
0039 

0036 

0037 

0038 

0039 
0060 
0061 
0062 
0063 
0066 

0065 

0066 
0067 
0066 


C 


c 

c 


SUBROUTINE STUFF ( J/X/SA/ VAX/ VR/Al/ OX/ OSR/ 0T1/K7/K8/ RONS/ TAH33/TF 

» 6/ TRM33/ TRM36/ NORAD/ N0AX/NZ/N6 ) , , . 

THIS SECTION CONFUTES THE COMPENSATORY TERNS USED IN EQSt Sib, 

,0, F^i?SlV R BY l Z£RS!!!!i!.. 0 F76 AND 0P81 IN STUFF., 
to FORMAT l ' 01 V BY ZERO 0PT6 ™ ' ) 

111 FORMAT! ' NEC SORT 0P76 AND 0P81 IN STUFF' 

112 FORMAT! 'NEC SORT. .>•• 0P76 AND 0F83 IN STUFF', 

TRM33«0# JTRM36«0, J TRM39«0t JTRN36»0, 

0P30»,l./(2.*3al6159,,*0Tl 

0P30A«DP30*0X 

OP306«0P30*DSR 


33 

36 


933 

917 

909 

910 


00 22 Kll»l/R8 

lF,Kll,EQ.1.0R.Kil.EQ.K8,0P71»0.5 

AKll«Kll-l 

Xl»OX*AKil 

SR3»0SR*AK11 

|F(Kll.CE,K8,Xl-Xl-0X*.5 
IP ,Kll.GE|K8 ,SR3»5R3«D5R*.9 

OP77»O.j0P78»O.|OP66«O.i0P83»O, 

DO 21 K9»l/K7 
AK9«K9-1 
T1«0T1*AK9 
I F(K0N3,EQ.1,RS»1 • 

IF (KUN3.EQ.2,RS«SR 

0P100«( (X-Al ,**2*RS**2-2a*RS*SR3*C0S,Tl,+SR3**2,6*3 

IF (K0N3 , OT, 1 ) CD TO 33 

IF(0P100,953/953/36 

IFIOPIOO, 917/909/36 

QP73*50RT , DP100 , 

IF(NURAD,EQ,2,G0 TO 1 

0F76»VAX,Rll!*(R5-SR3*C0StTl,I95R3/OP73 

CO TO 910 
0P76 *0 a i 00 TO 910 
0P76"0aJWRI1E, 6/311, JCO TO 910 
0P76-0. JWRlTE(6/309> 

OP77»OP77*OP76 


C 

C 


1 

911 


UNIT l,a|aaa«aOP«CAX/OXaa.aaaaaaaa 

I F ( RUNS , EQ . 1 ) 00 TO 938 
IF (NUAX, EQ . 2 ) 00 TO 32 
JFlOPlOO.lE.Oa ,00 TO 911 
0P81«VAX(KH)*,X-AL)6SR3/0P73 


60 TO 912 
PP81*0, 


TRM36 





0049 

912 

0090 

23 

0091 

e 

0092 

c 

0093 

938 

0096 


0099 


0036 

92 

0037 


0099 


0039 

96 

0060 


0061 


0062 


0063 

997 

0064 


0069 

932 

0066 

918 

0067 

913 

0068 

914 

0069 


0070 

C 

0071 


0072 


0073 


0074 

2 

0073 

6 

0076 

9 

0077 

C 

0078 


0079 


0080 

919 

0081 


0082 

80 

0083 


0084 


0083 

913 

0086 

916 

0087 

C 

0088 

21 

0089 

C 

0090 


0091 


0092 


0093 

3 

0094 


0093 


0096 


0097 

C 

0098 

22 

0099 

C 

0100 


0101 


0102 


0103 

4 

0104 


0103 


0106 


0107 

C 

0108 

931 

0109 



OP83»0P83*CIP8t 
l P (KQN5 , GT ■ 1 ) CO TO 32 

UNIT 3 TRM33 «, 

IPUl.GE.ADGO TO 932 

0P8Z»( < Al-Xl )**2*SR4*2-2.*SR*CQS(T1)U. )**3 
lP(0P82)932/932/36 

QP82* ( (X“Xl)**2*SR**2«2,*SR*COS(Tl )♦!« )**3 
IP<NOAX, £9,2)00 TO 2 
IP{ OP 82)918/913/36 
OP75«S9RT(OP82) 

IP (K0N3 , UE • 1 ) CO TO 93T 
0P76»VR(K11)*<X-X1)/0P73 
00 TU 914 

0P76«VR<Kll)*(AL-Xl)/0P73 
00 TO 914 
0P76»0i;00 TO 914 
0P76*0.jWRITE(6/312)jG0 TO 91* 
QP76«0,JWRITE<6/310) 

OP78«OP78*OP76 
IP (NORAOiEQa 2 )00 TO 21 

IPIK0N3.1T, 2)00 TO 919 
IP (0P82 .IE .0, )G0 TO 919 
GO TO 3 

IP<0P82)913/913/6 

0P73«S9RT(0P82) 

continue 

UN J T 4..,...,.0P-CR/DR 

OP83»VR<KU)*(SR-COS(T1))/DP73 
00 TU 916 

0PB2«{ (X-X1)**2-2.*C0S(T1)*2,>**3 
IP<0P82)915, 915/80 
OP73«S9RT(DP82 1 

0P83«VRIKU J*( 1,-COS (Tl))/OP75 

00 TO 916 

0PB3»0. 

0P84»0P84*0P83 

CONTINUE 

IP(NORAO,EQ,2)CO TO 3 
TRN3*"TRM34*OP7l*OP84 
TRM3S»TRM33*UP7l*0P77 
Jp (NQAX, EQ, 2 ) 00 TO 22 
TRM33«TRH33*0P71*0P78 
IP (K0N3.1T ,2)00 TO 22 
TRM36„TRM36*OP71*OPB3 

0P71«1. 

IP (NORAD, EQ ,2)00 TO 4 

TRM3A«TRH34*0P30A 

TRM33»TRM33*QP30B 

I P ( NUAX , EQ , 2 ) CO TO 931 

TRM33«TRM33*0P30A 

I F (KON 3 ,IT ,2)00 TO 931 

TRH36»TRM36*OP30B 






RETURN 

END 





CD 

0001 


i 

i 

SUBROUTINE COOP ( X> SR/ ICONS* A4/B4* A3* 63/ SUM2* SUM3* SUH4/SUM5* DTI* OP 

0002 

t 

*T1*N0RAD*N0AX* Il*Al*NZ»N4) 

000* 

C THIS SECTION COMPUTES THE FREESTREAM-lNlET INDUCED VELOCITY T6KH5* 

0004 

307 

FORMAT! 1 01V BY ZER0....0PI3 EO....IN COOPt) 

0005 


SUM2*0.jSUM3*0.j5UM4*0.;SUM3«0, 

0006 


DO 6 1*1/11 

0007 


AI* 1*1 

0006 


IPU|EQ,1,0R,:.EQ.I1)0P14*0.5 

0009 


T1*DT1*AI 

0010 


IF(K0N5. £0.2)00 TO 2 

0011 


A4*AL**2*SR**2 

0012 


64«>2.*SR*C0S(T1) 

0013 


A3»X**2*1 . 

0014 


B3--2.*C0StTl) 

0015 


00 TO 3 

0016 

2 

A4«X**2+SR**2 

0017 


64*»2«*SR*C0S(T1 ) 

0016 


AS*A4 

0019 


63*64 

0020 

3 

CONTINUE 

0021 


Al*A3jBl»83jA2*A4jB2*B4 

0022 


IF (NOAX. EQ. 2 )00 TO 1 

0023 

929 

AS*4,«A2*B2**>2 

0024 


IF ( A5.EQ.0i )C0 TO 933 .... 

0025 


0P13«(1./(4,*A2-B2**2) )*<4.*S0RT(A2)-U2.9B2+4.*A2)/S0RT(A24B24l 

0026 

t 

)) 

0027 


0P16*0P14«0P13«0T1 

0028 


GO TO 935 

0029 

933 

0P16*0.iWRlTE(6/3O7) 

0030 

935 

I F ( K0N5 , OE .2)00 TO 930 

0031 


5UM3*SUM3*0P16 

0032 


SRR-1. 

0033 


GO TO 931 

0034 

930 

SUM5*SUH5*0P16 

0035 

1 

I F ( NORAD i EQ i 2 ) CO TO 939 

0036 


SRR*SR 

0037 

931 

0P1*4,*A1-61**2 

0036 


lF(ABS(OPl)ilE..001 ) OP 1*. 0000000000000000000000000000 

0039 


1FI0P1.EQ.0, |GD TO 900 

0040 


0P9*I1./«4.*A1-61*«2)) . . ... ... ,, 

0041 


0P10»( t2,*Bl**2-4.*Al*2.*Al*6l)/ {SQRT (A1*B 1*1, ) ) )-2,*Bi*SQRT ( Al ) 

0042 


UP7*0P9*(4i*S0RT(Al)«((2.»Bl*4,*Al)/(SQRT(Al*Bl + l.) ) ) )*SRR 

0043 


OP20i«S0RT(Al)+81/2, 

0044 


IF(OP20l)900/900/200i 

0045 

2001 

0Pll*ALQ0((SQRT(Al + Bl*i.)*Bl/2.*l.)/0P201) 

0046 


0P6*C0S(T1 >*<OP9*OPIO+OP11) 

0047 


QP12*OP7“OP8 

0046 


OP15«OP14*QP12*OT1 

0049 


GO TU 903 

0050 

900 

0P15*O. 

0051 

903 

I F (KQN5.GE .2)00 TO 932 

0052 


SUM2*SUM2*0P15 

0053 


CO TO 939 

0054 

932 

SUM4»SUM4*0P15 

0055 

939 

CONTINUE 

0056 

6 

OP 14 * l • 

0057 

3003 

RETURN 

0056 


END 



APPENDIX II 


PROCEDURE FOR INPUTTING 
COMPUTER PROGRAM VARIABLES 
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The Program Options 

LAST - This quantity is inputted as 2 for normal operations and inputted 
as 1 if only the nacelle generation sections of the program are 
to be operated. 

NQTET - Inputted as 1 if reingestion is to be studied; otherwise 0. 

NOCROS - Inputted as 1 if cross -ingestion is to be studied, otherwise 0. 

NSPEED - Inputted as 0 if only one aircraft speed is to be studied. Other- 
wise the value of this quantity is the number of aircraft speeds 
to be studied . 

NOCARD - Inputted as 0 if nacelle generation sections are to be used. 

If this quantity is set equal to 1 , the singularity-produced 
velocities are read in on computer cards , rather than determined 
in the program. 

NOPEST - Inputted as 1 if the program is to determine whether or not the 
entraining portion of the reversed jet lies within the pre-entry 
streamtube. Otherwise, inputted as 0. 

Of the terms NOJET, NOCROS, and NOPEST, only one can be 

non-zero during a particular study. 

The Program Inputs 

EP - The degree of accuracy to which the no-flow condition is satisfied . 
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DELTAX - The dimensionless incremental value of x used in the 

streamline computational section. If this quantity is inputted 
as positive the streamlines will be calculated in the direction 
of the freestream flow. If "DELTAX” is negative, the opposite 
direction is used. 

The Dynamic Input 

U8A - The dimensional velocity of the aircraft. If more than one velocity 
is to be studied, the highest velocity is inputted here. 

VELRAT - The inlet-to-freestream velocity ratio. If more than one 

velocity is to be studied, the lowest velocity ratio is inputted 
here . 

VELJET - The dimensional reversed jet velocity. 

The Geometric Input 

ALA - The dimensional nacelle length. 

Circle - This quantity expands the concept of the Maximum 

Penetration Point from a point to an area. "CIRCLE" is 
inputted as the dimensionless diameter of a circle with the 
center at the Maximum Penetration Point. If this concept is 
not to be used, the quantity is inputted as zero. 
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The remaining terms in this section have the same meaning as in 
the main body of this report, but must be inputted in dimensional form. 
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